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PREFACE 


The  Central  Technical  Unit,  Hydrology  Branch,  Engineering  Division,  Soil 
Conservation  Service  (SCS),  U.S.  Department  of  Agriculture  (USDA) ,  uses 
electronic  data  processing  (EDP)  to  make  volume -duration -probability  (VDP) 
analyses  of  runoff  for  selected  watersheds.  Annual  maximum  streamflow 
data  for  periods  of  1,  3 ,  1 ,  15>  30,  60,  90,  120,  183,  and  27^  days  are 
furnished  to  SCS  by  the  U.S.  Geological  Survey  (USGS) . 

The  EDP  programs  described  in  this  technical  paper:  (l)  Transfer  USGS 
taped  data  directly  to  punched  cards  ready  for  processing,  (2)  perform 
all  the  statistical  computations  needed  to  obtain  the  maximum  annual  vol¬ 
ume  of  flow  at  each  selected  percent  chance  for  the  10  listed  periods, 
and  (3)  plot  the  computed  results  on  log  normal  probability  charts.  The 
principal  program  (0872)  uses  the  two -parameter  gamma  distribution  to 
compute  the  0-  to  99”Percent  chance  events  unless  the  gamma  shape  factor 
exceeds  51  j  then  it  uses  the  log  normal  distribution.  The  program  devel¬ 
ops  'a  wide  range  of  statistics  that  apply  to  most  cumulative  probability- 
distribution  analyses  of  natural  data  in  the  earth-science  field. 

The  programs  have  been  developed  jointly  by  the  Mathematical  and  Survey 
Applications  Section,  Washington  Data  Processing  Center,  Statistical 
Reporting  Service  (SRS),  and  the  Central  Technical  Unit,  Hydrology  Branch, 
Engineering  Division,  SCS,  USDA.  They  are  written  in  FORTRAN  II  for 
processing  through  the  IBM  707^>  IBM  1401,  and  EAI  Dataplotter  Model 
3300  equipment  available  in  the  Washington  Data  Processing  Center 

This  technical  paper  includes  source  programs,  tabular  data,  instructions 
for  preparing  input  data,  and  some  examples  of  application.  It  is  pri¬ 
marily  a  reference  to  the  analyses  being  made  in  the  Soil  Conservation 
Service.  It  is  intended  also  as  a  guide  for  making  similar  analyses  by 
other  technicians  who  are  familiar  with  electronic  data  processing. 
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HYDROLOGY  STUDY--A  MULTIPURPOSE  PROGRAM  FOR  SELECTED 
CUMULATIVE  PROBABILITY-DISTRIBUTION  ANALYSES 

By  W.  H.  Sammons,  hydraulic  engineer.  Soil  Conservation  Service 


IBM  7 07 A  PROGRAM  ABSTRACT 

Title:  Hydrology  study — a  multipurpose  program  for  selected  cumulative 

probability-distribution  analyses 

Program:  0872 

Job  No.:  920603-0010 

Date:  September  1963 

Purpose:  To  compute  various  statistics  from  hydrologic  or  other  data. 

Equipment  specifications: 

(a)  10,000  word  7070/72/7A  series 

(b)  Seven  tape  units: 

One  program  tape 

Two  input  tapes  (units  10  and  21) 

Four  output  tapes  (units  11,  12,  20,  and  22) 

(c)  Compiled  without  floating-point  hardware 

Source  language:  FORTRAN 

Accuracy:  Single  precision  floating  point 
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PROGRAM  0872  WRITEUP 

A.  Title:  Hydrology  study — a  multipurpose  program  for  selected  cumulative 
probability -distribution  analyses 

B.  Authors:  Elgin  G.  Fry,  Head,  Mathematical  and  Survey  Applications 
Section,  Washington  Data  Processing  Center,  SRS,  USDA,  and  William  H. 
Sammons,  Central  Technical  Unit,  SCS,  USDA 

C.  Description:  This  program  reads  specified  control  input  data  (described 
in  later  sections)  that  give  information  to  the  computer  concerning 

the  number  of  lines  plus  the  beginning-  and  ending-line  numbers  of  any 
desired  sets  of  100  lines  for  the  specified  input.  A  maximum  of  10 
sets  can  be  requested.  A  line  is  one  horizontal  line  of  input  data  on 
one  IBM  80 -column  card.  Only  79  columns  are  used;  see  data  card  format. 

A  maximum  of  100  lines  can  be  processed  for  each  station.  Each  line 
card  is  prepared  as  outlined  on  SRS -DP -102,  Multiple -Card  Layout  Form, 
page  11.  The  program  assumes  that  the  decimal  is  placed  just  before 
the  last  position  of  each  field.  Punching  of  decimals  in  the  data 
cards  overrides  the  program  and  the  data  field  will  be  read  incorrectly. 

Input  data  as  mean  discharge  for  10  selected  periods,  i.e.,  1,  3*  7> 

15  ,  30  ,  60  ,  90>  120,  183,  and  274  days,  are  arranged  in  descending  order 
of  magnitude  of  the  original  data  values .  The  values  are  converted  to 
log  to  base  10  and  then  to  inches .  The  following  values  are  written 
as  output  for  each  stream  gage: 

Original  values  from  high  to  low 

Log  to  base  10  of  original  values 

Mean  discharge  converted  to  inches  of  runoff 

Additional  values  required  for  automatic  data  plotting  are  outlined 
in  E. 

D.  Principal  statistics: 

1.  Number  of  lines  per  station  (N) 

where  X>0;  2  ^  N  ^  100 

N  is  printed. 

2.  Coefficient  of  variation  (CV) 

Used  with  arithmetic  normal,  extreme  value,  etc.,  for  small 
values  of  CV;  assumes  no  outliers  in  the  sample  size,  N. 
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2  N2  2  (X2)  -  N  (SO2 

CV  =  (N  -  l)  (H)y 

.  ' .  CV  =  V  CV2 


CV  is  printed. 

^  3.  Arithmetic  mean  (X  BAR) 

Used  with  arithmetic  normal,  extreme  value,  gamma,  etc.,  distri¬ 
butions;  assumes  no  outliers  in  the  sample  size,  N. 

Y  _  ^  X 

X  “  N 


X  BAR  is  printed. 

4.  Log  to  base  e  of  geometric  mean  (LN  G) 

Used  with  log  normal,  log  extreme  value,  gamma,  etc.,  distribu¬ 
tions;  outliers  may  or  may  not  be  present — in  general  none  should 
be  present. 


X  log  X. 
L°ge  G  =  — f-i 


LN  G  is  printed. 
j  5*  Geometric  mean  (G) 

Used  with  log  normal,  log  extreme  value,  gamma,  etc.,  distribu 
tions;  assumes  no  outliers  in  the  sample  size,  N. 

G  =  e^loge  G) 


G  is  printed. 

6.  Ratio  of  X  to  G 


Used  in  maximum-likelihood  solutions  of  log  normal,  log  extreme 
value,  gamma,  etc.,  distributions;  assumes  no  outliers  in  the 
sample  size,  N. 


X/G  = 


X  X 
N 


gT 


RATIO  X/G  is  printed  for  X/G. 


h 


7-  Coefficient  of  variation  (maximum  likelihood)  (CV*) 

Used  if  the  coefficient  of  variation  is  larger  than  l/3.  Used  with 
arithmetic  normal,  log  normal,  extreme  value,  log  extreme  value, 
gamma,  etc.,  distributions;  assumes  no  outliers  in  the  sample 
size,  N. 


CV*2  =  [  (X/G)2  -  1  ] 
. * .  CV*  =  V  CV*2 


CV*  is  printed. 

8.  Gamma -shape  parameter  (y) 

Used  in  two-  and  three -parameter  gamma  distributions  and  in  items 
9  and  10. 


Y  =  Loge(X/G) 

(a)  If  0  <  Y  <  0.5772 

y  =  Y  O.5OOO876  +  0.1648852Y  -  0.0544274Y2) 

(b)  If  0.5772  <  Y  <  17.0 

8 . 898919  +  9.059950Y  +  0.9775575Y2 

Y  Y(17 .79728  + 11.9^8477Y  +  Y^) 

Values  of  Y  and  GAMMA  are  printed.  If  Y  >  17-0,  the  remaining 
calculations  for  the  respective  set  are  not  made. 

9.  LN  S.D. 

Used  if  a  biased  estimate  of  the  log  standard  deviation  by  the 
maximum -likelihood  method  is  needed;  assumes  no  outliers  in  the 
sample  size,  N. 


LN  S.D.  =  y  2  Y 
LN  S.D.  is  printed. 

10.  LN  S.D.* 


Used  if  an  unbiased  estimate  of  the  log  standard  deviation  by 
the  maximum-likelihood  method  is  needed;  assumes  no  outliers  in 
the  sample  size,  N. 
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LN  S.D.*  =  7  ,  N  v 

MU-1' 


LN  S.D.*  is  printed. 

11.  LN  S.D.** 

Used  if  an  unbiased  estimate  of  the  log  standard  deviation  by 
the  method  of  moments  is  needed,  i.e.,  log  normal  distribution; 
assumes  no  outliers  in  the  sample  size,  N.  The  solution  is 
obtained  as  in  item  2. 

LN  S.D.**  =  CV  |LN  G| 

z 1  1 


where  Z  =  log  X. 

toe  i 

LN  S.D.**  is  printed,  but  the  definition  is  not  printed  on 
the  computer  print-out  sheets.  See  item  18. 

12.  Square  root  of  gamma  (SQL  GAMMA) 

Used  in  items  15,  1 6,  and  17 . 

SQR  GAMMA  =  V~y" 


SQL  GAMMA  is  printed. 

13.  Scale  parameter  (BETA) 
Used  in  item  15. 


B  =  (X/y) 

BETA  is  printed.  BETA  is  calculated  in  the  original  units .  If 
the  drainage  area  is  not  zero,  BETA  is  coded  as  BETA*,  inches. 

14.  Scale  parameter  converted  to  inches  (BETA*) 

B*  =  [  (D)  (0.037190083)  (B)  ]/M 

where  D  =  duration  in  days 

M  =  drainage  area  in  square  miles 

if  M  =  0,  B*  =  B 

B*  is  in  inches  of  depth  and  is  printed  as  BETA*. 


b 


15 .  BETA  times  SQR  GAMMA,  the  standard  deviation  in  the  original 
units 


B  *  S(GA)  =  b\T7 
B  *  S(GA)  is  printed. 

16.  BETA*  times  SQR  GAMMA,  the  biased  standard  deviation  in  inches  of 
depth 


B*  (S(GA)  =  (B *)V~7" 

B*  (S(GA)  is  printed. 

17.  U(P,IT)  times  BETA*  times  SQR  GAMMA 

U(B*(S(GA)  =  (U)  (B*)  (V^T) 

where  U=  Value  of  f  (r  ,P) 

P=  (y  -  1) 

If  P  is  between  two  levels  of  the  shape  parameter  (y  —  l),  the  P 
value  for  selected  probabilities  of  I’  (greater  than  I'  is  a 
selected  probability  level)  is  determined  by  linear  interpolation. 

See  tables  1  and  2  for  I’  =  0.0;  0.2;  1;  2;  4;  10;  20;  50;  80;  90; 
and  99  percent  (greater  than) . 

Table  1  is  used  if  P  =  —0.95(0.05)4  (pp.  12-17). 

Table  2  is  used  if  P  =  —0.05(0. 5)74(l)l64  (pp.  18-30). 

Table  3  supplements  table  1  (pp.  31-4l)« 

The  computer  program  uses  an  abbreviated  table  that  includes  all 
of  table  1  and  values  for  P  =  4.5(0.5)38(1)50  from  table  2.  If 
P  >  50  (or  GAMMA  >  51);  the  log  normal  is  computed.  See  items 
11  and  18. 

U(B*(S(GA)  is  printed  for  the  11  computed  points  listed. 

18.  Log-normal  solution  (calculated  only  if  GAMMA  >  51*0) 

X  =  [  (e(log  G  +  K  S  X)  )  (D)  (0.37190083)  ]/M 

fee 

N2L(Z.2)  -  N(EZ. )2 
<CVz>2  =  (N  -  l)  (ZZi)2 

where  M  =  drainage  area  in  square  miles 


D  =  duration  in  days 
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=  reduced  variate  for  selected  probability  levels  (item  19) 

Z.  =  log  (X.)  is  the  log  to  base  e  of  the  transformed  variate 
1  e  1 

S  =  log  standard  deviation  or  LN  S.D.** 

(CV  )2  is  the  unbiased  coefficient  of  variation  of  the  log  trans¬ 
formed  variate  and  is  not  printed.  See  item  11. 

19.  Other  cumulative  probability  distributions 

Reduced  variates  for  selected  distributions 
for  given  return  periods 


Return 

period 

^Proba¬ 

bility 

y 

value 

K 

0* 

valQe 

K 

value 

1.01 

99 

-1.52718 

-1.6408 

-2.5263 

1.05 

95 

-1.09719 

-1*3055 

-1.6449 

1.11 

90 

-0.83^03 

-1.1003 

-1.2816 

1*25 

80 

-O.47588 

-0.8211 

-0.8416 

2 

50 

O.3665I 

-0.1643 

0.0000 

5 

20 

1.49994 

0.7194 

0.8416 

10 

10 

2.25037 

1.3045 

1.2816 

25 

4 

3.19853 

2.0438 

1.7507 

50 

2 

3  *  90194 

2.5923 

2.0537 

100 

1 

4.60015 

3.1367 

2.3263 

500 

0.2 

6.21561 

4.3947 

2.8782 

105 

0.001 

11.51292 

8.5265 

4.2649* 

PM 

0.00 

21.87455 

15.82576 

5.250** 

6.015 

*  Not  used  in  this  program. 

Probability  (percentage  greater  than) . 

**  Seven  digits  +5  considered  for  probability;  this  is  a  graphic 
relationship. 


Other  sources  give  these  reduced  variate  values  different  mag¬ 
nitudes.  For  seven  9's  and  a  5,  5*0400  to  5*3270,  5*250  is  used 
as  a  compromise.  For  nine  9's  and  a  5,  5*485  to  6.110,  6.015  is 
used  as  a  compromise  since  from  a  graphic  plotting  5*250  and  6.015 
are  consistent. 


The  return  period  (T)  is  defined  as 


where  I’  =  the  probability  (greater  than) 

For  the  extreme  value  theory ,  y  value 

y  =  —  in  (—  In  I) 
where  I  =  the  probability  (less  than) 
fn  =  the  natural  logarithm 

The  equation  for  X.  computed  for  selected  y  values  (y  is  a  function 
of  the  return  period) 


=  u  +  (l/a)  y 


where  u  =  the  mode 

(l/a)  =  the  logarithmic  rate  of  increase  (slope) 

y  =  the  reduced  variate  for  selected  probabilities 

or  for  the  log  extreme  value  distribution 

X.  =  log  _1  [  +  (l/a)ig  y  ] 

or  in  terms  of  K  values 
g 


X.  =  X  +  K  S 
1  -  g  x 

where  X  =  the  arithmetic  mean,  item  3 

K  =  the  Gumbel  reduced  variate  for  extreme  value  distribution 
g 

S  =  the  standard  deviation  in  the  original  units,  which  can 
be  calculated  from  item  2  or  7>  depending  on  the  size 
of  the  coefficient  of  variation 

S  =  CV  X  or  CV*  X 


x 

and  in  terms  of  the  log  extreme  value  distribution 

X.  =  log  _1  [  X„  +  K  S  ,  ] 
l  &  fg  -  g  xjfcg 

where  X  .  =  item  4  converted  to  log  base  10 

fg 


xig 


=  item  9;  10,  or  11  converted  to  log  base  10 
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K  value  is  used  with  the  normal  and  log  normal  distributions, 
n  ° 

For  the  normal 


X.  =  X  +  K  S 
1  —  n  x 


For  the  log  normal 

X.  =  log""1  (X.  +  K  S  .  ) 

1  °  v  ig  -  n  xJL g' 

Both  K  and  K  assume  population  limitations  where  N  -  oo. 
g  n 


Sample  distributions  of  K  and  K  can  be  considered  if  needed. 

g  n 


E.  Preparation  of  input  data:  Data  must  be  punched  on  80-column  cards 
(only  79  columns  are  used  in  this  program)  in  the  exact  layout  speci¬ 
fied  on  the  multiple -card  layout  form.  Columns  1  and  2  are  used  to 
identify  each  line  of  data  arranged  from  largest  to  smallest  (YR  on 
the  print-out  exhibit,  pp.  91 “93) •  Identification  can  range  from  00 
to  99.  There  are  10  fields  or  columns  of  input  data — 7  fields  of  8 
digits  and  3  fields  of  7  digits.  Column  80  was  not  used  in  this 
version.  It  was  reserved  for  a  binary  code  for  the  IBM-7090  computer 
program,  which  preceded  this  IBM-707^  program. 


F.  Header  card  or  parameter  card:  Data  must  be  punched  on  80-column 
cards  in  the  exact  layout  specified  on  the  multiple -card  layout  form. 
Columns  1  through  10  are  for  station  identification.  Columns  11 
through  20  are  for  drainage  area  (square  miles).  The  period  of  record 
is  always  equal  to  the  total  or  set  1  and  must  be  coded  as  001  in 
columns  21  through  23  to  total  N  in  columns  2k  through  26.  Sets  2 
through  10  can  be  any  sequence  of  numbers  from  001  to  N  or  parts 
thereof.  It  is  possible  to  use  all  columns  from  21  through  80  if  10 
sets  of  computations  are  desired.  If  more  than  10  sets  are  needed,  a 
new  station  is  designated. 


G.  End-of-file  card:  All  9's  must  be  punched  in  columns  1  through  30  on 
an  80-column  card  as  shown  in  the  layout  on  the  multiple-card  layout 
form. 


H.  Data  plotter  preparation:  During  the  processing  of  each  set  or  subset, 
the  values  of  CAFX(i)  (conversion  of  mean  discharge  to  inches)  are 
converted  to  log^Q  and  written  on  tape  unit  20  to  be  used  as  input  to 
program  0911,  which  prepares  the  values  for  use  with  the  automatic 
data  plotter.  If  GAMMA  ^  51*0,  the  values  of  UB*  V  y  are  converted 
to  log, n  and  written  immediately  after  the  CAFX(i)  values.  If  GAMMA 
is  >  51.0,  the  log  normal  solutions  for  X.  are  converted  to  log  10  and 
written  out. 
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The  ordinate  ranges  from  a  lower  limit  of  10  4  =  0.0  or  0.0001  units  to 
no  set  upper  limit.  The  abscissa  is  the  normal  reduced  variate  (Kn) 
where  a  lower  limit  of  —4.0000  =  0.0  and  +6.0000  =  0.9999  or  1.0.  Plot¬ 
ting  of  the  observed  ordered  data  is  based  on  normal  order  statistics 
for  the  abscissa  (Kn)  values.  They  can  be  plotted  in  line  or  dot  mode. 
The  11  computed  points  are  given  in  items  17,  18  and  19  of  D.  The 
ordinate  is  the  sum  of  the  log  to  base  10  of  the  observed  or  computed 
data  and  +4.000000,  which  is  then  multiplied  by  0.1.  A  scale  range  must 
be  selected  manually  for  the  data  plotted  from  the  ADP  print -out.  See 
programs  0910  and  0911  in  the  exhibits . 

Print -outs  of  tables  1,  2,  and  3  and  the  different  programs  follow. 


MULTIPLE-CARD  LAYOUT  FORM 
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Table  3 •  — Percentage  points  of  the  chi-square  distribution-^ 

v  =  0. 1(0.1)  10.0 
y  =  0.05(0.05)  5-0 

p  =  -0. 95(0. 05)  4.0 


The  values  in  the  body  of  Harter’s  table  2  are  the  values  of  chi  square 
that  correspond  to  these  values  of  v  (degrees  of  freedom,  reading  across, 
then  down)  for  P  =  .0001,  .0005,  .0010,  .0050,  .0100,  .0250,  .0500,  .1000 
(,1000)>  .9000,  .9750,  .9900,  .9950,  .9990,  .9995,  and  .9999  (probability 
(less  than)  and  I’  =  100 (l  —  P)).  Table  1,  Percentage  points  of  incom¬ 
plete  gamma -function  ratio,  I'  (U,P)  was  based  on  this  table.  Table  3  (10  pp.) 
supplements  table  1  for  probabilities  not  listed  in  table  1.  Basic 
formulas  for  this  table  follow: 


X2/ 2  =  Xj,/3  =  uVy" 
v  =  2y  3  =  X/y 
Example  1:  Given  y  =  0.25 


v  =  2y  =  0.50 


=  ugVV" 

X.  =  (X./p)(p) 

P  =  Y  -  1  =  -0-75 
x2  =  0.875289E  01 


I'  =  0.1  or  P  =  0.9990  u  =  8.75289 


Example  2:  Given  y  =  0*50  v  =  2y  =  1.0  p=— 0.50 

I’  =  50.0  or  P  =  0,5000 

X2/2  =  uV0.50  x2  =  0.454936 

0.494236  =  ^0.56 

0.227468  =  u(  0.70  7IO678H) 
u  =  0.32168833  or 

0.321689  as  read  dirfectly  from  table  1, 
page  1  of  5 


3/  After  Harter,  H.  L.  More  tables  of  the  incomplete  gamma -function  ratio 
and  of  percentage  points  of  the  chi-square  distribution.  U.S.  Air  Force 
Aerospace  Res.  Lab.. Tech.  Rpt.  64-123.  92  pp.  1964.  (For  sale  Office 

oT  Technical  Services,  U.S.  Dept.  Commerce,  $2.25.) 
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percentage  POINTS  OF 

THE  CHI-SQUARE 

DISTRIBUTION 

p  V 

0.1 

0.2 

0.3 

0.4 

.000  1 

0. 1 1 6893E— 79 

0.121 46 1 E— 39 

0.271457E— 26 

0. 13051 OE-19 

0. 

.0005 

0.11 i 478E— 65 

0. 1 1861 4E-32 

0.  124023E-21 

0.407843E— 16 

0. 

.0010 

0.1 16893E-59 

0. 121461E-29 

0.125999E— 19 

0. 13051 OE-14 

0. 

.0050 

0.11 1 478E— 45 

0.1 18614E-2 2 

0.575663E-1 5 

0.407843E— 1 1 

0. 

.0100 

0.1 16893E-39 

0. 121461 E— 1 9 

0 . 584837E— 1 3 

0. 13051 0E-09 

o. 

.0250 

0.106313E-31 

0. 1 1 5834E— 1 5 

0 . 263007E- 1 0 

0. 12745 1 E— 07 

0. 

.0500 

0.11 i 478E— 25 

0. 1 18614E-12 

0.2671 99E— 08 

0 • 407843E— 06 

0. 

•  1  000 

0.1 16893E-19 

0. 121461 E— 0  9 

0.271457E— 06 

0. 13051 0E-04 

o. 

.2000 

G.122571E-13 

0. 124376E— 06 

0.275757E-04 

0.41 770  3E— 03 

o. 

.3000 

0.4  07579E— 1 0 

0.71 721 7E-05 

0.411 706E-03 

0.31 7558E-02 

o. 

.4000 

0. 1 28525E— 07 

0. 1 27368E-03 

0.280520E-02 

0. 134391E-01 

o. 

.5000 

0.11 1 478E— 05 

0. 1 18678E— 02 

0. 124696E-01 

0.414927E— 01 

0. 

.6000 

0.427387E-04 

Q.73689GE— 02 

0.425976E— 01 

0. 1 0602 1 E— 00 

0. 

.7000 

0.9331 27E— 03 

0 . 348556E— 0 1 

0. 123216E-00 

0.242075E— 00 

0. 

.8000 

0 • 1 35640E— 0 1 

0.1 38780E-00 

0.326574E-00 

Q.527087E  00 

0. 

.9000 

0. 152634E-00 

0.532309E  00 

0 . 889597E  Ou 

0.  120980E  01 

0* 

.9500 

0 . 53 1 865E  00 

0. 1 16087E  01 

0.1 65 1 1 7E  01 

0.206105E  01 

o. 

.9750 

0. 1 13435E  01 

0.1 9558 1 E  01 

0 . 254506E  01 

0.302221E  01 

0. 

.9900 

0.217525E  01 

0.31 7696E  01 

0.3861 19E  01 

0.440461E  01 

0. 

.9950 

0 • 308925E  01 

0.41 89 1 2E  01 

0.492727E  01 

0 . 550939E  01 

o. 

.9990 

0.547292E  01 

0.672735E  01 

0 .75552 1 E  01 

0.820393E  01 

0. 

.9995 

0 . 657533E  01 

0 . 787676E  01 

0 . S73335E  01 

0.940409E  01 

0. 

.9999 

0.924820E  01 

0.1Q6325E  02 

0.1 15418E  02 

0.122539E  02 

0* 

P  V 

0.6 

0.7 

0.8 

0.9 

.000  1 

0 .647288E- 1 3 

0.536368E-1 1 

0. 148307E-09 

0. 197223E— 08 

0* 

.0005 

0.138356E-10 

0.532746E— 09 

0.829062E-08 

0. 705058E-07 

0. 

.0010 

0. 139454E-09 

0.386016E-08 

0.468988E-07 

0 . 328988E— 06 

0. 

.0050 

0 • 2 980 7 9E— 07 

0.38341 OE-06 

0.262173E-05 

0.11 761  IE-04 

o. 

.0100 

G.300445E-06 

0.27781 1 E— 05 

0.1 4S308E— 04 

0 . 548796E— 04 

0. 

.0250 

0.6371 36E— 05 

0.380826E-04 

0. 146566E-03 

0.42051 0E-03 

0. 

.0500 

0 .642207E-04 

0.275963E-03 

0.829307E-03 

0. 196319E-02 

0. 

.  1  000 

0 • 64  744 9E— 03 

0.200084E— 02 

0 . 469775E  — 02 

0.91 833 1 E— 02 

0. 

.2000 

0 .654068E-02 

0. 145652E— 01 

0.267845E-01 

0.433564E— 01 

o. 

.3000 

0 • 254533E— 0 1 

0 • 469482E— 0 1 

0.750838E— 01 

0. 109174E-00 

0* 

.4000 

0.674796E— 01 

0. 109267E-00 

0. 1 58724E-Q0 

0.21 4367E— 00 

0. 

.5000 

0. 146262E-00 

0.21 4739E— 00 

0.290 1 56E-QO 

0.37C657E— 00 

o. 

.6000 

0 .282505E-00 

0.383777E— 00 

0.489505E-0O 

0 . 596032E  00 

o. 

.7000 

0 • 5 1 3 1 30E  00 

0.653661E  00 

0.794514E  00 

0.934828E  00 

0. 

.8000 

0.920148E  00 

0. 1 10836E  01 

0.1 29 1 1 4E  01 

0.146896E  01 

o. 

.9000 

0.176962E  01 

0.202139E  01 

0 • 225969E  01 

0 . 2487 1 OE  01 

0« 

.9500 

0 • 274470E  01 

0 . 304392E  01 

0 • 332392E  01 

0.35888  IE  01 

o. 

.9750 

0.380053E  01 

0.413667E  01 

0 . 444943E  01 

0 . 474398E  01 

0* 

.9900 

0.527882E  01 

0.565317E  01 

0.600019E  01 

0 . 632602E  01 

0. 

.9950 

Q.644114E  01 

0 . 683889E  01 

0.720'709E  01 

0.75524  IE  01 

0. 

.9990 

0 • 923787E  01 

0.967813E  01 

0.1 00852E  02 

0.1  04666E  02 

o. 

.9995 

0.104728E  02 

0.109278E  02 

0.1 13485E  02 

0.117427E  02 

0. 

.9999 

0.1 33893E  02 

0.1 3873 1 E  02 

0.1 43205E  02 

0 • 1 47399E  02 

0. 

1  of  10 

0.5 

134994E-15 
843712E-13 
134994E-1 1 
843712E-09 
134994E— 07 
527320E-06 
Q43715E-05 
1 3500 1 E— 03 
216177E-02 
1 09826E— 0 1 
350448E-01 
873476E— 0 1 
1 884 1 2E— 00 
374696E— 00 
7261 70E  00 
150078E  01 
242023E 
343324E 
486777E 
600362E 
875289E 
997154E 
128566E 
1  .0 

157080E-07 
392699E— 06 
1 57080E— 05 
392704E-04 
157088E-03 
982069E-03 
3932 1 4E— 02 
157908E-01 
641848E-01 
1 48472E— 00 
274996E— 00 
454936E— 00 
708326E  00 
107419E 
164237E 
270554E 
384146E 
502389E 
66349QE 
787944E 
108276E 
121 157E 
151367E 


01 

01 

01 

01 

01 

01 

02 


01 

01 

01 

01 

01 

01 

01 

02 

02 

02 


33 


2  of  10 

PERCENTAGE  POINTS  OF  THE  CHI-SQUARE  DISTRIBUTION 


V 

1  •  1 

1  .2 

1 .3 

1  .4 

1  .5 

0001 

0 

•  86 1 549E  — 07 

0 

•  357 1 63E— 0  6 

0 

. 1 19364E- 

05 

0. 

33674 7E 

-05 

0 

. 8295 1 5E 

-05 

0005 

0 

•  1 60744E-05 

0 

.5221 77E-05 

0 

. 1 41 976E- 

04 

0. 

33561  IE 

-04 

0 

• 709238E 

-04 

0010 

0 

• 5 668 4 3E— 05 

0 

.  165781E-04 

0 

•412422E- 

04 

0. 

90341 3E 

-04 

0 

• 178722E 

-03 

0050 

0 

•  1 05762E-03 

0 

•242391 E-03 

0 

.49061 3E- 

03 

0. 

900572E 

-03 

0 

• 152864E 

-02 

0100 

0 

•372989E-03 

0 

•  76967 OE  — 03 

0 

. 1 42556E- 

02 

0. 

242525E 

-02 

0 

• 385450E 

-02 

0250 

0 

.  1 97446E-02 

0 

•354744E-02 

0 

. 504496E— 

02 

0. 

899664E 

-02 

0 

•1311 30E 

-01 

0500 

0 

•697388E— 02 

0 

•  1 12897E— 0 1 

0 

. 1 70363E- 

01 

0. 

243264E 

-01 

0 

. 332328E 

-01 

1000 

0 

•247335E— 0 1 

0 

.361 209E-0 1 

0 

•499625E- 

01 

0. 

66291 OE 

-01 

0 

. 8498 1 OE 

-01 

2000 

0 

.890340E-0 1 

0 

. 1 1 7607E-00 

0 

. 1 49589E- 

00 

0. 

184677E 

-00 

0 

•222586E 

-00 

3000 

0 

. 1 92276E-00 

0 

•23997BE-00 

0 

.291 046E- 

00 

0. 

345040E 

-00 

0 

•401589E 

-00 

4000 

0 

•  33 966 7E—  00 

0 

•407645E— 00 

0 

•  478359E- 

00 

0. 

551360E 

00 

0 

.  626292E 

00 

5000 

0 

•542073E  00 

0 

•631404E  00 

0 

.722448 E 

00 

0. 

814847E 

00 

0 

. 908334E 

00 

6000 

0 

•81 973 1 E  00 

0 

•  93 181  BE  00 

0 

.  1 04430E 

01 

0. 

1 15700E 

01 

0 

• 126978E 

01 

7000 

0 

•  1 2 i 244E  01 

0 

•134950E  01 

0 

.  1 48538E 

01 

0. 

162013E 

01 

0 

•  1 75379E 

01 

8000 

0 

•  1 8 i 1 90E  01 

0 

• 197798E  01 

0 

.2141 OOE 

01 

0. 

230129E 

01 

0 

•245912E 

01 

9000 

0 

• 29 i 647E  01 

0 

•312101E  01 

0 

• 332003E 

01 

0. 

351426E 

01 

0 

•370426E 

01 

9500 

0 

•408400E  01 

0 

•431602E  01 

0 

•  4544  78E 

01 

0. 

476524E 

01 

0 

.  498020E 

01 

9750 

0 

•529176E  01 

0 

• 554955E  01 

0 

. 579674E 

01 

0. 

60405 1 E 

01 

0 

.62758 IE 

01 

9900 

0 

•692986E  01 

0 

•721318E  01 

0 

•748660E 

01 

0. 

775148E 

01 

0 

•800890E 

01 

9950 

0 

•  81914  7E  01 

0 

• 849097E  01 

0 

•677960E 

01 

0. 

905943E 

01 

0 

•933102E 

01 

9990 

0 

.  1  1 1717E  02 

0 

•  1 1501 7E  02 

0 

. 1 16198E 

02 

0. 

121275E 

02 

0 

•  1 24262E 

02 

9995 

0 

•124712E  02 

0 

•128122E  02 

0 

•  1 31 406E 

02 

0. 

134587E 

02 

0 

. 137672E 

02 

9999 

0 

• 155151E  02 

0 

. 158780E  02 

0 

. 1 62277E 

02 

0. 

165660E 

02 

0 

•  168944E 

02 

V 

1  *6 

1 . 7 

1 .8 

1  .9 

2.0 

000  1 

0 

•  1 82997E-04 

0 

•U66620E-04 

0 

. 688306E— 

04 

0. 

120566E 

-03 

0 

.20001 OE 

-03 

0005 

0 

. 1 36826E— 03 

0 

•444862E-03 

0 

.41 1560E- 

03 

0. 

65621  IE 

-03 

0 

. 1 00025E 

-02 

0010 

0 

•  32 54 4 7E— 03 

0 

• 553490E— 03 

0 

.8891 73E- 

03 

0. 

136143E 

-02 

0 

. 200 1 OOE 

-02 

0050 

0 

•243471E-02 

0 

•367955E— 02 

0 

•532263E- 

02 

0. 

74204  IE 

-02 

0 

. 100251E 

-01 

0100 

0 

.57961 5E-02 

0 

• 83270 7E—0 2 

0 

•  1  15163E- 

01 

0. 

1 54238E 

-01 

0 

•201007E 

-01 

0250 

0 

.  1&2B39E-01 

0 

• 245789E— 0 1 

0 

• 320466E- 

01 

0. 

407274E 

-01 

0 

• 506356E 

-01 

0500 

0 

•  4379tDlE-01 

0 

• 560263E— 0 1 

0 

.6991 69E- 

01 

0. 

8545Q3E 

-01 

0 

.  102587E 

-00 

1  000 

u 

. 1 05964E-00 

0 

.  1291  37E-r00 

0 

•154393E- 

00 

0. 

181624E 

-00 

0 

•210721E 

-00 

2000 

0 

•263057E-00 

0 

•305858E-00 

0 

.  350783E- 

00 

0. 

397647E 

-00 

0 

. 446287E 

-00 

3000 

0 

• 460 3 8 4E— 00 

0 

•  52 1 1 6 1 E  00 

0 

. 583702E 

OO 

0. 

647818E 

00 

0 

•713350E 

00 

4000 

0 

•702«71E  00 

0 

•  78066  7E  00 

0 

. 860094E 

00 

0. 

940398E 

00 

0 

.  1021ET5E 

01 

5000 

0 

•100270E  01 

0 

• 109779E  01 

0 

.  1  19349E 

01 

0. 

128968E 

01 

0 

•  1 38629E 

01 

6000 

0 

• 138254E  01 

0 

•  1 49524E  01 

0 

.  160783E 

01 

0. 

172028E 

01 

0 

• 183258E 

01 

7000 

0 

.  188643E  0  1 

0 

.20181  IE  01 

0 

.2 14668E 

01 

0. 

22788  IE 

01 

0 

.240795E 

01 

8u00 

0 

•  26 l 4  74E  01 

0 

•276836E  01 

0 

•2920 1 5E 

01 

0. 

307028E 

01 

0 

•321888E 

01 

9000 

0 

•  3890i>2E  01 

0 

•407343E  01 

0 

•  425332E 

01 

0. 

443049E 

01 

0 

.460517E 

01 

9500 

0 

•519029E  01 

0 

•539605E  01 

0 

•559792E 

01 

0. 

579629E 

01 

0 

•599146E 

01 

9750 

0 

•650S39E  01 

0 

• 672989E  01 

0 

• 694983E 

01 

0. 

716566E 

01 

0 

•  737776E 

01 

9900 

0 

•825976E  01 

0 

•830477E  01 

0 

. 874454E 

01 

0. 

697959E 

01 

0 

.921 034E 

01 

9950 

0 

• 959554E  01 

0 

• 985375E  01 

0 

.  1 01063E 

02 

0. 

103538E 

02 

0 

. 105966E 

02 

9990 

0 

•  1271  fc>9E  02 

0 

• 130005E  02 

0 

.  1 32778E 

02 

0. 

135492E 

02 

0 

• 138155E 

02 

9995 

0 

•140675E  02 

0 

. 143604E  02 

0 

.  1 46466E 

02 

0. 

149269E 

02 

0 

•152018E 

02 

9999 

0 

. 172139E  02 

0 

•175255E  02 

0 

.  1 7830 1 E 

02 

0. 

181283E 

02 

0 

. 1842Q7E 

02 

3b 


3  of  10 

PERCENTAGE  POINTS  OF  THE  CHI-SQUARE  DISTRIBUTION 


V 

2. 1 

2.2 

2.3 

2.4 

2.5 

GOO  1 

0 

•316675E- 

03 

0 

• 48 1 566E- 

03 

0 

•  7  07  046E- 

03 

0 

. 1 00666E- 

02 

0 

.  139492E- 

02 

0005 

0 

. 1 46697E— 

02 

0 

•  206089E- 

02 

0 

.266725E- 

02 

0 

.385157E- 

02 

0 

.50591 7E- 

02 

00  1  0 

0 

•  283963E  — 

02 

0 

•390934E- 

02 

0 

.524 1 69E- 

02 

0 

. 666743E— 

02 

0 

.  881 588E  — 

02 

0050 

0 

. 1 3 i 838E  — 

01 

0 

. 1 69385E- 

0  1 

0 

•213252E- 

01 

0 

. 263750E— 

01 

0 

.321 134E- 

01 

0100 

0 

•255868E- 

01 

0 

.31921 6E- 

0  1 

0 

. 39 1 249E- 

01 

0 

.4721 74E- 

01 

0 

. 562 121 E- 

01 

0250 

0 

.61 761 7E- 

01 

0 

•741 649E- 

01 

'0 

•677765E- 

01 

0 

.  1  02602E- 

00 

0 

. 1 18623E- 

00 

0500 

0 

•  1 2 i 286E  — 

00 

0 

.141 504E- 

00 

0 

. 163188E- 

UU 

0 

. 186290E- 

00 

0 

.21 0756E- 

00 

1  000 

0 

.24 i S62E- 

00 

0 

.2741 OGE- 

00 

0 

.30821 OE- 

ou 

0 

. 343797E- 

00 

0 

•380789E- 

00 

2000 

0 

• 496559E- 

00 

0 

.548335E 

00 

0 

.601 500E 

ou 

0 

. 655952E 

00 

0 

.71 1 599E 

00 

3000 

0 

.780 1 60E 

00 

0 

. 848 1 3 1 E 

00 

0 

•917160E 

00 

0 

•987156E 

00 

0 

e 1 05804E 

01 

4000 

o 

. 1 1 0375E 

01 

0 

. 1 18660E 

01 

0 

.  127012E 

01 

0 

• 135425E 

01 

0 

.  143894E 

01 

5000 

0 

. 1 48327E 

01 

0 

. 1 58055E 

01 

0 

. 167610E 

01 

0 

. 1 77567E 

01 

0 

. 187385E 

01 

60  0  0 

0 

. 1 9447 1 E 

0  1 

0 

. 205667E 

0  1 

0 

.2 16846E 

01 

0 

.226007E 

01 

0 

•239150E 

01 

7000 

0 

• 253634E 

01 

0 

. 266403E 

0  1 

0 

.2791 07E 

01 

0 

.291 749E 

01 

0 

. 3Q4333E 

01 

6000 

0 

•  33660fc>E 

01 

0 

.3511 95E 

01 

0 

• 365662E 

01 

0 

. 38001 8E 

01 

0 

. 394268E 

01 

9000 

u 

•477756E 

01 

0 

.49479  IE 

0.1 

0 

.511 630E 

01 

0 

• 528292E 

01 

0 

. 544788E 

01 

9500 

0 

•61 6373E 

01 

0 

.63733 3E 

0  1 

0 

• 656047E 

01 

0 

• 674533E 

01 

0 

.69280BE 

01 

9750 

0 

. 75864  5E 

01 

0 

. 779202E 

01 

0 

.799471 E 

01 

0 

.619474E 

01 

0 

. 839230E 

01 

9900 

0 

. 9437 1 8E 

01 

0 

. 966042E 

0  1 

0 

. 988035E 

01 

0 

•  1 00972E 

02 

0 

.1031 12E 

02 

9950 

0 

. 1 08352E 

02 

0 

.  1  1 Q700E 

02 

0 

.  1  1301  IE 

02 

0 

. 1 15290E 

02 

0 

•  1 17538E 

02 

9990 

0 

. 140770E 

02 

0 

. 143341E 

02 

0 

. 1 4587 1 E 

02 

0 

. 148364E 

02 

0 

. 150822E 

02 

9995 

0 

. 1 5471 7E 

02 

0 

. 157370E 

02 

0 

. 159981E 

02 

0 

. 162553E 

02 

0 

.  165089E 

02 

9999 

0 

. 187077E 

02 

0 

. 1 89899E 

02 

0 

. 1 92675E 

02 

0 

. 195409E 

02 

0 

.  198104E 

02 

V 

2.6 

2.7 

2.8 

2.9 

3.0 

O0o  1 

0 

.  188710E- 

02 

0 

.24990 1E- 

02 

0 

•324678E- 

02 

0 

•414666E- 

02 

0 

•521483E- 

02 

0005 

0 

•  65 1 48 1 E  — 

02 

0 

•824242E- 

02 

0 

. 1 02648E- 

01 

0 

. 1 26Q36E- 

01 

0 

. 1 52790E- 

01 

0010 

0 

.111  1 4  7E- 

01 

0 

•  1 37896E- 

0  1 

0 

. 168644E- 

01 

0 

. 20360 5E- 

01 

0 

•242976E- 

01 

0050 

0 

.38561 5E- 

01 

0 

. 457355E— 

0  1 

0 

. 536462E- 

01 

0 

• 623083E— 

01 

0 

.717218E- 

01 

0100 

0 

.661 181E- 

01 

0 

. 769326E- 

0  1 

0 

.  886606E- 

01 

0 

. 1Q1296E- 

00 

0 

. 1 14832E- 

00 

0250 

0 

. 1 3581 6E- 

00 

0 

.  1 54157E- 

00 

0 

. 1 7361 9E- 

oo 

0 

. 1941 75E- 

00 

0 

.21 5795E- 

00 

0500 

0 

•23654 1 E— 

00 

0 

•263589E- 

00 

0 

.291 653E- 

00 

0 

.321287E- 

00 

0 

. 351 846E— 

00 

1  000 

0 

.4191 12E- 

00 

0 

.  456695E— 

00 

0 

. 4994  7  3E  — 

06 

0 

. 541385E 

00 

0 

. 584374E 

00 

2000 

0 

•  76836 1 E 

00 

0 

•826162E 

00 

0 

. 884938E 

oo 

0 

•944626E 

00 

0 

. 10051 7E 

01 

3000 

0 

. 1 12975E 

01 

0 

. 1 2022 1 E 

0  1 

0 

.  127539E 

01 

0 

•  1 3492 1 E 

01 

0 

. 142365E 

01 

4000 

0 

. 15241 3E 

01 

0 

. 1 6097  ^E 

0  1 

0 

. 1 69567E 

01 

0 

. 1 76234E 

01 

0 

. 186917E 

01 

5000 

0 

. 1 97200E 

01 

0 

• 207030E 

0  1 

0 

•216674E 

ol 

0 

• 226730E 

01 

0 

• 236597E 

01 

6000 

0 

•250277E 

01 

0 

.261 386E 

01 

0 

•272479E 

01 

0 

. 283556E 

01 

0 

. 2946 1 7E 

01 

7000 

G 

.31 6862E 

01 

0 

. 329340E 

01 

0 

.341 768E 

01 

0 

.354150E 

01 

0 

. 366487E 

01 

6UG0 

0 

•408420E 

01 

0 

• 42246 1 E 

01 

0 

. 436455E 

01 

0 

. 450347E 

01 

0 

•464168E 

01 

9000 

o 

•56 i 1 30E 

01 

0 

. 57732  7E 

0  1 

0 

• 593389E 

01 

0 

. 609324E 

01 

0 

.625139E 

01 

9500 

0 

.71 0886E 

01 

0 

. 728780E 

0  1 

0 

*  74650  3E 

01 

0 

. 7640 64 E 

01 

0 

•781473E 

01 

9750 

0 

• 858755E 

01 

0 

•  67606  7E 

0  1 

0 

.8971 76E 

01 

0 

.916097E 

01 

0 

. 934840E 

01 

9900 

0 

•  1 05226E 

02 

0 

•  1 0731 5E 

02 

0 

. 1 0936 1 E 

06 

0 

.  1  1  1425E 

02 

0 

. 1 13449E 

02 

9950 

0 

. 1 19757E 

02 

0 

. 121949E 

02 

0 

•  1  24 1  1 6E 

02 

0 

.  126260 E 

02 

0 

• 128382E 

02 

9990 

0 

. 1 53247E 

02 

0 

. 1 55642E 

02 

0 

.  158009E 

02 

0 

• 160348E 

02 

0 

. 162662E 

02 

9995 

0 

. 1 6759 1 E 

02 

0 

. 1 7006 1 E 

02 

0 

.  1 72502E 

02 

0 

. 174914E 

02 

0 

.177300E 

02 

9999 

0 

•20Q763E 

02 

0 

•203367E 

02 

0 

•  205979E 

02 

0 

• 208542E 

02 

0 

.21 1075E 

02 

35 


4  of  10 

PERCENTAGE  POINTS  OF  THE  CHI-SQUARE  DISTRIBUTION 


p  V 

3.  1 

3.2 

3.3 

3.4 

3.5 

•  0001 

0 

•646729E— 

02 

0 

•791963E 

-02 

0 

•958702E- 

02 

0 

.  1  1 4840E- 

01 

0 

• 136246E- 

01 

•  0005 

0 

.  183096E- 

01 

0 

.217126E- 

-0 1 

0 

•255038E- 

01 

0 

• 296973E— 

01 

0 

. 343058E- 

01 

.0010 

0 

• 286929E— 

01 

0 

.33561  7E- 

-01 

0 

•3e91 75E- 

01 

0 

.44771 7E- 

01 

0 

.51 1 340E- 

01 

.0050 

0 

•81 89 1 5E- 

01 

0 

.9281 81 E 

-0  1 

0 

. 1 04500E- 

00 

0 

•  1 1 6934E- 

00 

0 

.1301 16E- 

00 

.0100 

0 

. 129289E- 

00 

0 

. 1 44568E 

-00 

0 

. 1 60745E- 

OO 

0 

. 1 77777E- 

00 

0 

•  1 95651 E- 

00 

.0250 

0 

• 238452E- 

00 

0 

.2621 1 6E 

-00 

0 

•286760E- 

OO 

0 

•31 2355E- 

00 

0 

•338873E- 

00 

.0500 

0 

•3834o7E- 

00 

0 

•41616SE 

-00 

0 

•449850E- 

00 

0 

. 484495E— 

00 

0 

•520069E 

00 

.  1000 

0 

• 6283o9E 

00 

0 

•67338 1 E 

00 

0 

•71 9304E 

00 

0 

•  7661 16E 

00 

0 

•81 3778E 

00 

.2000 

u 

. 1 06653E 

01 

0 

. 1 12865E 

0  1 

0 

.  1  19149E 

01 

0 

. 125502E 

01 

0 

.131920E 

01 

.3000 

6 

.  1 498C7E 

0  1 

0 

. 157422E 

0 1 

0 

.  165027E 

01 

0 

. 17268  IE 

01 

0 

.  180379E 

01 

.4000 

0 

. 1 95634E 

01 

0 

•204382E 

01 

0 

.213160E 

01 

0 

•221964E 

01 

0 

•230795E 

01 

.5000 

0 

•246474E 

0  1 

0 

• 256359E 

01 

0 

•266252E 

01 

0 

.276152E 

01 

0 

•286059E 

01 

.6000 

0 

•  305662E 

01 

0 

•31 6694E 

01 

0 

.32771 OE 

01 

0 

.33671 3E 

01 

0 

.349703E 

01 

.7000 

0 

. 378763E 

01 

0 

. 39 1 038E 

0  1 

0 

•403256E 

01 

0 

•415437E 

01 

0 

•427563E 

01 

.8000 

0 

•477906E 

0  1 

0 

.491 580E 

0 1 

0 

.505 1 88E 

01 

0 

•51 8735E 

01 

0 

• 532222E 

01 

.9000 

0 

•640S41E 

01 

0 

. 656437E 

01 

0 

.671 931E 

01 

0 

. 687329E 

01 

0 

. 702636E 

01 

.9500 

0 

• 798738E 

01 

0 

.81 5868E 

01 

0 

• 832869E 

01 

0 

.649749E 

01 

0 

• 8665 1 2E 

01 

.9750 

0 

.95341 6E 

01 

0 

.971832E 

0  1 

0 

.  990099E 

01 

0 

. 1 00822E 

02 

0 

•  1  0262 1 E 

02 

.9900 

0 

•  1 1 5453E 

02 

0 

. 1 1 7439E 

02 

0 

.  1  19407E 

02 

0 

. 121359E 

02 

0 

. 123296E 

02 

.9950 

o 

. 1 30482E 

02 

0 

• 132563E 

02 

0 

. 1 34625E 

02 

0 

. 136669E 

02 

0 

•  138696E 

02 

.9990 

o 

. 1 64VD3E 

02 

0 

. 167220E 

02 

0 

.  169466E 

02 

0 

• 171692E 

02 

0 

. 173899E 

02 

.9995 

0 

. 1 7966  IE 

02 

0 

. 181998E 

02 

0 

. 184313E 

02 

0 

. 18660 7E 

02 

0 

• 188680E 

02 

.9999 

0 

•21 3582E 

02 

0 

.21 6063E 

02 

0 

•218520E 

0* 

0 

. 220954E 

02 

0 

. 22336SE 

02 

p  V 

3.6 

3.7 

3.8 

3.9 

4.0 

•  OOO  1 

0 

. 160220E- 

01 

0 

•  1 86887E 

-0  1 

0 

.21 6365E- 

01 

0 

•248764E- 

01 

0 

•284185E- 

■01 

.0005 

0 

•393405E- 

01 

0 

.4481 1 5E 

-01 

0 

• 507273E- 

01 

0 

. 570954E- 

01 

0 

.639220 E- 

01 

.0010 

0 

.500 128E- 

01 

0 

•654145E 

-01 

0 

• 733444E- 

01 

0 

.81 8067E- 

01 

0 

. 908040E- 

•01 

.0050 

0 

•  1 44039E- 

00 

0 

. 1 58696E 

-00 

0 

. 1 74079E- 

00 

0 

. 190180E- 

00 

0 

•206989E- 

■00 

.0100 

0 

•21 435 1 E- 

00 

0 

. 233863E 

-00 

0 

.254 1 70E- 

00 

0 

. 275258E- 

00 

0 

.2971 09E- 

•00 

.0250 

0 

.366289E- 

00 

0 

. 394577E 

-00 

0 

.42371 OE- 

00 

0 

. 453665E- 

00 

0 

. 48441 9E- 

00 

.0500 

0 

•556536E 

00 

0 

. 593865E 

00 

0 

• 632025E 

00 

0 

. 670987E 

00 

0 

•71 0723E 

00 

•  1  000 

0 

•862253E 

00 

0 

♦911 506E 

00 

0 

.961 505E 

ou 

0 

•  101222E 

01 

0 

. 106362E 

01 

.2000 

0 

• 138399E 

0  1 

0 

. 1 44937E 

01 

0 

. 151 53 1 E 

01 

0 

. 158179E 

01 

0 

• 164678E 

01 

.3000 

0 

. 188 12  IE 

01 

0 

. 1 95902E 

01 

0 

•203722E 

01 

0 

.21 1579E 

01 

0 

•219470E 

01 

.4000 

0 

• 239650E 

01 

0 

•248527E 

0 1 

0 

•257426E 

01 

0 

. 266346E 

01 

0 

•275284E 

01 

.5000 

0 

• 29597 1 E 

01 

0 

•305889E 

01 

0 

.31581  IE 

01 

0 

. 325738E 

01 

0 

•335669E 

01 

•  60U0 

u 

• 360679E 

01 

0 

•371643E 

01 

0 

•  382595E 

01 

0 

• 393534E 

01 

0 

• 404463E 

01 

.7000 

0 

•  439696E 

0  1 

0 

.451 777E 

01 

0 

. 463828E 

01 

0 

• 475850E 

01 

0 

•  487843E 

01 

.8000 

0 

• 545653E 

01 

0 

. 559030E 

0  1 

0 

• 572356E 

01 

0 

. 585632E 

01 

0 

•598862E 

01 

.9000 

0 

.71 7856E 

01 

0 

. 732994E 

0  1 

0 

. 748052E 

01 

0 

. 763034E 

01 

0 

. 777944E 

01 

.9500 

0 

.8831 65E 

01 

0 

.89971 4E 

0  1 

0 

•916161E 

01 

0 

•  93251 3E 

01 

0 

•948773E 

01 

.9750 

0 

. 1 04407E 

02 

0 

. 1 061 81E 

02 

0 

.  107942E 

02 

0 

•  1 09693E 

02 

0 

.111 433E 

02 

.9900 

0 

.  12521 7E 

02 

0 

.  127124E 

02 

0 

•129018E 

02 

0 

. 130899E 

02 

0 

.  132767E 

02 

.9950 

0 

.  1 40707E 

02 

0 

. 142702E 

02 

0 

. 1 44683E 

02 

0 

• 146650E 

02 

0 

. 148603E 

02 

.9990 

0 

•  1 76086E 

02 

0 

. 1 78257E 

02 

0 

•  1 804 1 OE 

02 

0 

• 182547E 

02 

0 

• 184668E 

02 

.9995 

0 

.19i 134E 

02 

0 

. 1 9337  OE 

02 

0 

•  1 95588E 

02 

0 

. 197789E 

02 

0 

.  199974E 

02 

.9999 

0 

•  225756E 

02 

0 

•228127E 

02 

0 

•230479E 

02 

0 

. 2328 1 2E 

02 

0 

•235127E 

02 

36 


percentage  points  of 


v 

4. 1 

4.2 

000  1 

0 

•322722E- 

0  1 

0 

•364463E- 

0  1 

0005 

0 

.712126E- 

01 

0 

•769712E- 

01 

0010 

0 

. 1 00336E- 

00 

0 

•  1 1 041 OE- 

00 

0050 

0 

•224496E- 

00 

0 

*  24269 1 E- 

00 

0100 

0 

.31971 OE- 

00 

0 

•343044E- 

00 

0250 

0 

•51 6947E 

00 

0 

•546228E 

00 

0500 

0 

•  75 i 2C7E 

00 

0 

.79241 5E 

00 

1000 

0 

. 1 1 i 569E 

0  1 

0 

. 1 1 6839E 

0  1 

2000 

0 

. 1 7 i 625E 

01 

0 

. 1 78420E 

01 

3000 

0 

•227394E 

01 

0 

. 235350E 

0  1 

4000 

0 

• 28424 1 E 

0  1 

0 

. 2932 1 6E 

01 

5000 

0 

• 345604E 

01 

0 

.  355542E 

01 

6000 

0 

.41 5360E 

01 

0 

. 426286E 

0  1 

7000 

0 

.49981 OE 

01 

0 

.511 750E 

0  1 

6000 

0 

*61 2046E 

01 

0 

•625186E 

0  1 

9000 

0 

.7927o5E 

01 

0 

•807559E 

0  1 

9500 

0 

. 964945E 

01 

0 

•981032E 

0  1 

9750 

0 

*  1  1 3 1 62E 

02 

0 

•  1  1 4882E 

02 

9900 

0 

. 134623E 

02 

0 

.  1 36468E 

02 

9950 

0 

. 1 50B43E 

02 

0 

. 1 52470E 

02 

9990 

0 

. 1 86775E 

02 

0 

•  1 88867E 

02 

9995 

0 

•  202 1 43E 

02 

0 

. 204297E 

02 

9999 

0 

• 237426E 

02 

0 

. 239709E 

02 

V 

4.6 

4.7 

000  1 

0 

• 564938E- 

01 

0 

.623744E- 

0  1 

0005 

0 

. 1 14742E- 

00 

0 

. 1 24876E- 

00 

00  1  0 

0 

. 1 56067E- 

00 

0 

. 1 688 1 5E- 

•00 

0050 

0 

•  322 1  A  1 E- 

00 

0 

.34361 7E- 

•00 

0100 

u 

,443403E- 

00 

0 

•470175E- 

•00 

0250 

0 

• 664475E 

00 

0 

. 720220E 

00 

0500 

0 

»  964  0 1 5E 

00 

0 

. 1 008506 

0  1 

1000 

0 

. 1 36509E 

01 

0 

. 144064E 

0  1 

2000 

0 

.206029E 

0  1 

0 

.21 3031 E 

0  1 

3000 

u 

. 267465E 

01 

0 

. 275560E 

0  1 

4000 

0 

• 329272E 

0  1 

0 

. 338322E 

01 

5000 

0 

. 395325E 

01 

0 

»  405277E 

0  1 

6000 

0 

.46961  IE 

01 

0 

• 480668E 

0  l 

7000 

0 

•  55927 1 E 

01 

0 

.571 095E 

0  1 

6000 

0 

*  677346E 

0  1 

0 

•690292E 

01 

9000 

0 

•666036E 

U1 

0 

. U6U5 1 5E 

01 

950  0 

0 

. 1 04460E 

02 

0 

. 1 060316 

02 

9750 

0 

. 1 2l669E 

02 

0 

. 123345E 

02 

9900 

0 

•  1 4374 IE 

02 

0 

. 1 45535E 

02 

9950 

0 

. 160066E 

02 

0 

. 161936E 

02 

9990 

0 

. 1 971 04E 

02 

0 

. 1 991 32E 

02 

9995 

0 

.212775E 

02 

0 

.21 4862E 

02 

9999 

0 

• 248667E 

02 

0 

. 250897E 

02 

5  Of  10 

CHI-SQUARE  DISTRIBUTION 


4.3 

4.4 

4.5 

409487E- 

01 

0 

457865E- 

01 

0 

•509662E- 

01 

8720 13E- 

01 

0 

959054E- 

01 

0 

.  105085E- 

00 

121 020E- 

00 

0 

132168E- 

00 

0 

• 143850E- 

00 

261 564E- 

00 

0 

281 104E- 

00 

0 

• 30 1 300E- 

00 

367095E- 

00 

0 

391849E- 

00 

0 

•41 729QE- 

00 

581242E 

00 

0 

614967E 

00 

0 

. 649385E 

00 

83432 1 E 

00 

0 

876904E 

00 

0 

.920143E 

00 

1221 70E 

01 

0 

127560E 

01 

0 

• 133007E 

01 

1 85259E 

01 

0 

1 92 1 4 1 E 

01 

0 

• 199065E 

01 

243337 E 

01 

0 

251352E 

01 

0 

•259395E 

01 

302207E 

01 

0 

31 1213E 

01 

0 

•320235E 

01 

365484E 

01 

0 

375428E 

01 

0 

•385375E 

01 

437182E 

01 

0 

448068E 

01 

0 

•458944E 

01 

523666E 

01 

0 

535557E 

01 

0 

. 547425E 

01 

636285E 

01 

0 

651343E 

01 

0 

•664363E 

01 

822268E 

01 

0 

8369 17E 

01 

0 

•651506E 

01 

997038E 

01 

0 

101297E 

02 

0 

•  102882E 

02 

1 16592E 

02 

0 

1 18293E 

02 

0 

•  1 19985E 

02 

1 38302E 

02 

0 

140125E 

02 

0 

• 141938E 

02 

1 54386E 

02 

0 

156290E 

02 

0 

• 158183E 

02 

190946E 

02 

0 

19301  IE 

02 

0 

. 195064E 

02 

206437E 

02 

0 

208563E 

02 

0 

•21 0675E 

02 

241975E 

02 

0 

244227E 

02 

0 

. 246464E 

02 

4.8 

4.9 

1 

5.0 

686126E- 

01 

0 

752124E- 

01 

0 

•821774E-01 

135487E- 

00 

0 

146575E- 

00 

0 

• 15813SE- 

■00 

1 82091 E- 

00 

0 

1 95891 E- 

00 

0 

•210213E- 

00 

3657 1 6E- 

00 

0 

366428E- 

00 

0 

.41 1742E- 

•00 

49759QE- 

QO 

0 

525636E 

00 

0 

. 554298E 

00 

756602E 

Ou 

0 

793605E 

00 

0 

•831212E 

00 

1 05359E 

01 

0 

109925E 

01 

0 

•  1 14548E 

01 

14967  IE 

01 

0 

155327E 

01 

0 

•  1 61 03 1 E 

01 

220070E 

01 

0 

227144E 

01 

0 

•234253E 

01 

283680E 

01 

0 

291824E 

01 

0 

•299991 E 

01 

347385E 

01 

0 

35646 IE 

01 

0 

. 365550E 

01 

4 1 523 1 E 

01 

0 

425188E 

01 

0 

.435146E 

01 

491 51 6E 

01 

0 

502356E 

01 

0 

•513187E 

01 

582897E 

01 

0 

594680E 

01 

0 

•606443E 

01 

703204E 

01 

0 

716062E 

01 

0 

. 728928E 

01 

894940E 

01 

0 

909312E 

01 

0 

. 923636E 

01 

1 07596E 

02 

0 

109153E 

02 

0 

•  1 1 0705E 

02 

125012E 

02 

0 

126672E 

02 

0 

. 128325E 

02 

1 47320E 

02 

0 

149096E 

02 

0 

• 150863E 

02 

163800E 

02 

0 

165653E 

02 

0 

. 167496E 

02 

201 149E 

02 

0 

203155E 

02 

0 

.205150E 

02 

2 16937E 

02 

0 

2 1 9001 E 

02 

0 

.221 053E 

02 

253093E 

02 

0 

255277E 

02 

0 

•257448E 

02 

THE 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

0. 

0. 

0. 

0. 

0. 

0. 

0® 

0. 

o. 

0. 

0® 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

o. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


37 


6  of  10 

PERCENTAGE  POINTS  OF  THE  CHI-SQUARE  DISTRIBUTION 


V 

5.  1 

5.2 

5.3 

5.4 

5.5 

0001 

0 

.8951 05E- 

01 

0 

.972142E- 

0  1 

0 

•  105291E- 

00 

0 

1 13741E- 

00 

0 

• 122568E- 

00 

0005 

0 

. 1 701 74E- 

00 

0 

. 182683E- 

•00 

0 

. 195661 E- 

00 

0 

209106E- 

00 

0 

• 2230 1 5E- 

00 

0010 

0 

• 225051 E- 

00 

0 

. 240402E- 

00 

0 

• 256260E- 

00 

0 

272622E- 

00 

0 

•289481 E- 

00 

0050 

0 

•435647E- 

00 

0 

•460133E- 

00 

0 

•485190E- 

00 

0 

51 0807E 

00 

0 

•536975E 

00 

0100 

0 

*  583564E 

00 

0 

.613420E 

00 

0 

. 643853E 

00 

0 

674853E 

00 

0 

•  7064Q7E 

00 

0250 

0 

•869407E 

00 

0 

•908176E 

00 

0 

. 947504E 

00 

0 

987377E 

00 

0 

.  102778E 

01 

0500 

0 

. 1 19225E 

01 

0 

. 123955E 

01 

0 

. 128736E 

01 

0 

133568E 

01 

0 

•  1 38448E 

01 

1000 

0 

•  166781E 

01 

0 

• 172577E 

01 

0 

•  1784 16E 

01 

0 

184298E 

01 

0 

•  190221E 

01 

2000 

0 

• 24 1 396E 

01 

0 

. 248570E 

0  1 

0 

• 255776E 

01 

0 

263012E 

01 

0 

•270277E 

01 

3000 

0 

•308180E 

01 

0 

•31 6389E 

01 

0 

. 324620E 

01 

0 

332870E 

01 

0 

.341 1 40E 

01 

4000 

0 

• 374650E 

01 

0 

. 383762E 

0  1 

0 

. J92884E 

01 

0 

402016E 

01 

0 

•411 1 60E 

01 

5000 

o 

.4451 06E 

01 

0 

. 455068E 

01 

0 

•46503 1 E 

01 

0 

474996E 

01 

0 

• 484962E 

01 

6000 

0 

• 524009E 

01 

0 

•534824E 

01 

0 

.54563  IE 

01 

0 

556430E 

01 

0 

•567222E 

01 

7000 

0 

•618187E 

01 

0 

•629912E 

0  1 

0 

•641620E 

01 

0 

65331 OE 

01 

0 

. 664984E 

01 

8000 

0 

•  74 1 742E 

01 

0 

. 754526E 

01 

0 

. 767280E 

01 

0 

780006E 

01 

0 

•792704E 

01 

9000 

0 

.93791  IE 

0  1 

0 

•952140E 

0 1 

0 

. 966324E 

01 

0 

980465E 

01 

0 

•994563E 

01 

9500 

0 

•  1 12250E 

02 

0 

. 1 1 3790E 

02 

0 

. 1 15324E 

02 

0 

1 16852E 

02 

0 

•  1 18376E 

02 

9750 

0 

. 129971E 

02 

0 

. 131609E 

02 

0 

. 13324  IE 

02 

0 

134867E 

02 

0 

•  136486E 

02 

9900 

0 

• 152622E 

02 

0 

. 154373E 

02 

0 

. 1561 16E 

02 

0 

1 5785 1 E 

02 

0 

• 159580E 

02 

9950 

0 

. 169330E 

02 

0 

.171 1 56E 

02 

0 

. 1 72973E 

02 

0 

174782E 

02 

0 

• 176583E 

02 

9990 

0 

.207135E 

02 

0 

.2091 09E 

02 

0 

.21 1 074E 

02 

0 

213030E 

02 

0 

.214976E 

02 

9995 

0 

•223095E 

02 

0 

•225126E 

02 

0 

.'227146E 

02 

0 

229157E 

02 

0 

.231 158E 

02 

9999 

0 

•259608E 

02 

0 

•261756E 

02 

0 

•263893E 

02 

0 

26601 9E 

02 

0 

.2681 34E 

02 

V 

5.6 

5.7 

5.8 

5.9 

6.0 

000  1 

0 

•  13i 770E- 

00 

0 

.141 350E- 

00 

0 

•151 307E- 

00 

0 

161641 E- 

00 

0 

.  172352E- 

00 

0005 

0 

• 237385E- 

00 

0 

• 2522 1 3E- 

00 

0 

• 267495E- 

00 

0 

283228E- 

00 

0 

•29940BE- 

•00 

00  1  0 

0 

. 306834E- 

00 

0 

•324674E- 

•00 

0 

. 342997E- 

00 

0 

361 796E- 

00 

0 

. 381 067E- 

•00 

0050 

0 

. 563662E 

00 

0 

. 59092 1 E 

00 

0 

•61 868 1 E 

00 

0 

646952E 

00 

0 

•675727E 

00 

0100 

0 

• 738503E 

00 

0 

.771 1 30E 

00 

b 

• 804277E 

00 

0 

837934E 

00 

0 

. 872Q90E 

00 

0250 

0 

. 1 06870E 

01 

0 

.  1  1  1  0 1 3E 

0  1 

0 

•  1 15206E 

01 

0 

1  1 9447E 

01 

0 

• 123734E 

01 

0500 

0 

• 143376E 

01 

0 

• 148350E 

01 

0 

.  153370E 

01 

0 

158433E 

01 

0 

• 163538E 

01 

1  000 

0 

•  1 96184E 

01 

0 

.202185E 

01 

0 

.208225E 

ol 

0 

21 43Q 1 E 

01 

0 

•220413E 

01 

2000 

0 

•  27757 OE 

01 

0 

. 284890E 

0  1 

0 

•  292238E 

01 

0 

29961  IE 

01 

0 

•307009E 

01 

3000 

0 

•349428E 

01 

0 

. 357734E 

01 

0 

.  366058E 

01 

0 

374398E 

01 

0 

. 382755E 

01 

4000 

0 

.42031 3E 

01 

0 

. 429475E 

0  1 

0 

•  438646E 

01 

0 

447826E 

01 

0 

•457015E 

01 

5000 

0 

•  494930E 

01 

0 

•504899E 

0  1 

0 

•514869E 

01 

0 

524840E 

01 

0 

•534812E 

01 

6000 

o 

•576006E 

01 

0 

•588784E 

0  1 

0 

. 599555E 

01 

0 

61 031 8E 

01 

0 

.621 Q76E 

01 

7000 

0 

. 67664 1 E 

01 

0 

• 688282E 

0  1 

0 

.699907E 

01 

0 

71  151  BE 

01 

0 

.7231 14E 

01 

8000 

0 

.  805375E 

01 

0 

•81 8020E 

0  1 

0 

•  830640E 

01 

0 

843235E 

01 

0 

•855806E 

01 

9000 

0 

. 1 00862E 

02 

0 

. 1 02264E 

02 

0 

.  1 03662E 

02 

0 

105056E 

02 

0 

• 106446E 

02 

9500 

0 

. 1 19893E 

02 

0 

• 121406E 

02 

0 

.  1 229 1 4E 

02 

0 

12441 7E 

02 

0 

.  125916E 

02 

9750 

0 

. 138099E 

02 

0 

. 1 39706E 

02 

0 

. 141 308E 

02 

0 

1 42903E 

02 

0 

. 144494E 

02 

9900 

0 

• 161301E 

02 

0 

. 163015E 

02 

0 

. 164723E 

02 

0 

166424E 

02 

0 

. 1 68 1 1 9E 

02 

9950 

0 

•  1 78376E 

02 

0 

. 180162E 

02 

0 

. 1 81 940E 

02 

0 

18371 1 E 

02 

0 

. 185476E 

02 

9990 

0 

.216913E 

02 

0 

•218842E 

02 

0 

.220762E 

02 

0 

222674E 

02 

0 

• 224577E 

02 

9995 

0 

•233150E 

02 

0 

.2351 32E 

02 

0 

.2371 06E 

02 

0 

239071 E 

02 

0 

.241 028E 

02 

9999 

0 

• 27q240E 

02 

0 

•272335E 

02 

0 

•274420E 

02 

0 

276497E 

02 

0 

.278563E 

02 

38 


7  of  10 

percentage  points  of  the  chi-square  distribution 


P  V 

6.  1 

6.2 

6*3 

6.4 

6.5 

.000  1 

0 

. 1 83439E- 

00 

0 

. 1 94902E- 

00 

0 

.206740E- 

00 

0 

21 8950E- 

00 

0 

231 533E- 

00 

.0005 

0 

•31 6030E— 

00 

0 

. 33309 1 E- 

00 

0 

•350587E- 

00 

0 

368514E- 

00 

0 

386866E- 

00 

.0010 

0 

•400B04E- 

00 

0 

. 42 1001 E- 

00 

0 

.441 653E- 

00 

0 

462755E- 

00 

0 

484300E- 

00 

.0050 

0 

.  704995E 

00 

0 

•734749E 

00 

0 

• 764979E 

00 

0 

795677E 

00 

0 

826835E 

00 

.0100 

0 

.  906736E 

00 

0 

.94 186  IE 

00 

0 

. 977456E 

00 

0 

1 01 35  IE 

01 

0 

1 05002E 

01 

.  0250 

0 

.  128068E 

01 

0 

•  1 32447E 

0  1 

0 

. 1 36870E 

01 

0 

141 336E 

01 

0 

145845E 

01 

.0500 

0 

. 1 68686E 

01 

0 

.  1 73874E 

01 

0 

. 1 791 02E 

01 

0 

184369E 

01 

0 

189673E 

01 

•  1  000 

0 

•  226560E 

01 

0 

. 23274 1 E 

01 

0 

• 238954E 

01 

0 

245200E 

01 

0 

251478E 

01 

.2000 

0 

•31 443 1 E 

01 

0 

.32 1 878E 

0  1 

0 

.  32934  7E 

01 

0 

336839E 

01 

0 

344353E 

01 

.3000 

0 

.  39 i  1 28E 

01 

0 

.39951 6E 

01 

0 

.  407920E 

01 

0 

4 1 6338E 

01 

0 

424770E 

01 

.4000 

0 

•4662 1 3E 

01 

0 

.4754 1 8E 

01 

0 

.  48463 1 E 

01 

0 

493851 E 

01 

0 

503079E 

01 

.5000 

0 

• 544785E 

01 

0 

• 554759E 

0  1 

0 

•564734E 

01 

0 

57471 OE 

01 

0 

584687E 

01 

.6000 

0 

•  63 i 827E 

01 

0 

•  64257 1 E 

0 1 

0 

• 6533 1 OE 

01 

0 

664043E 

01 

0 

674769E 

01 

.7000 

0 

• 734695E 

01 

0 

. 746263E 

01 

0 

.75781 7E 

01 

0 

769358E 

01 

0 

780886E 

01 

.8000 

0 

•868354E 

01 

0 

•880878E 

0  1 

0 

•89338 1 E 

01 

0 

905862E 

01 

0 

918323E 

01 

.9000 

0 

•  1 07833E 

02 

0 

. 1 0921 7E 

02 

0 

. 1 10597E 

02 

0 

1 1 1 974E 

02 

0 

1 1 3347E 

02 

.9500 

0 

.  12741 OE 

02 

0 

. 128900E 

02 

0 

. 1 30385E 

02 

0 

131866E 

02 

0 

1 33343E 

0’2 

.9750 

0 

.  146079E 

02 

0 

. 1 47659E 

02 

0 

. 1 49234E 

02 

0 

150804E 

02 

0 

152369E 

02 

.9900 

0 

. 1 69808E 

02 

0 

. 1 71490E 

02 

0 

. 1 73167E 

02 

0 

174838E 

02 

0 

1 76504E 

02 

.9950 

0 

•  187234E 

02 

0 

•  1 66985E 

02 

0 

•  1 90730E 

02 

0 

192468E 

02 

0 

1 9420 1 E 

02 

.9990 

0 

•226474E 

02 

0 

.22 8362E 

02 

0 

•  23024  3E 

02 

0 

2321 1 7E 

02 

0 

233984E 

02 

.9995 

0 

•242977E 

02 

0 

.24491 7E 

02 

0 

.  246850E 

02 

0 

248775E 

02 

0 

250693E 

02 

.9999 

0 

.280621E 

02 

0 

. 28267 1 E 

02 

0 

.28471  IE 

02 

0 

286744E 

02 

0 

288768E 

02 

P  '  V 

6.6 

6.7 

6.8 

6.9 

7.0 

.000  1 

0 

.  244466E— 

00 

0 

. 257806E- 

00 

0 

.271 496E- 

00 

0 

285550E- 

00 

0 

299967E- 

00 

.0005 

0 

•  4  0564 1 E- 

00 

0 

.424834E- 

00 

0 

. 444440E- 

00 

0 

464455E- 

00 

0 

484875E- 

00 

.0010 

0 

•  506284E 

00 

0 

•526700E 

00 

0 

.551 545E 

00 

0 

57481 1 E 

00 

0 

598494E 

00 

.0050 

0 

. 858446E 

00 

0 

•  89050 1 E 

00 

0 

. 922992E 

00 

0 

95591 3E 

00 

0 

989256E 

00 

.0100 

0 

.  1 08697E 

01 

0 

. 1 1 2436E 

01 

0 

. 1 1 621 7E 

01 

0 

120040E 

01 

0 

1 23904E 

01 

.0250 

0 

. 1 50394E 

01 

0 

. 1 54984E 

0  1 

0 

. 1 5961 3E 

01 

0 

164281E 

01 

0 

168987E 

01 

.0500 

0 

•  1 9501 5E 

01 

0 

. 200393E 

0  1 

0 

• 205806E 

01 

0 

21 1 253E 

0  1 

0 

216735E 

01 

.  1  000 

0 

•  25  7766E 

01 

0 

.264124E 

0  1 

0 

•27049 1 E 

01 

0 

276887E 

01 

0 

28331  IE 

01 

.2000 

0 

• 35 i 8o8E 

01 

0 

. 359444E 

0  1 

0 

. 367021 E 

01 

0 

37461 7E 

01 

0 

382232E 

01 

.3000 

0 

•433216E 

0  1 

0 

.44 1 676E 

0 1 

0 

•450149E 

01 

0 

456635E 

01 

0 

467133E 

01 

.4000 

0 

•51231 5E 

01 

0 

•521557E 

0  1 

0 

. 530806E 

01 

0 

540061 E 

01 

0 

549324E 

01 

.5000 

0 

. 594664E 

01 

0 

•604643E 

0  1 

0 

•61 462 1 E 

01 

0 

624601 E 

01 

0 

63458 1 E 

01 

.6000 

0 

• 68549 1 E 

01 

0 

.696206E 

0  1 

0 

. 7069 1 6E 

01 

0 

71 762  IE 

01 

0 

72832 1 E 

01 

.7000 

0 

. 79240 1 E 

01 

0 

• 803904E 

0  1 

0 

•815396E 

01 

0 

826875E 

01 

0 

838343E 

01 

.8000 

0 

•930762E 

0  1 

0 

.9431 82E 

01 

0 

• 955582E 

01 

0 

967963E 

01 

0 

980325E 

01 

.9000 

0 

. 1 14718E 

02 

0 

. 1 1 6085E 

02 

0 

.  1  1 7450E 

02 

0 

1 1 881 2E 

02 

0 

1201 70E 

02 

.9500 

0 

•  1 3481 6E 

02 

0 

. 1 36286E 

02 

0 

. 1 3775 1 E 

02 

0 

13921 3E 

02 

0 

1  4067 1 E 

02 

.  9750 

u 

. 1 53930E 

02 

0 

. 1 55466E 

02 

0 

. 1 57037E 

02 

0 

1 56585E 

02 

0 

1601 28E 

02 

.9900 

0 

•1781 64E 

02 

0 

.  1  798 1 9E 

02 

0 

. 1 81 469E 

02 

0 

1831 1 3E 

02 

0 

1 84753E 

02 

.9950 

0 

. 195927E 

02 

0 

. 1 97648E 

02 

0 

.  1 99363E 

02 

0 

201 073E 

02 

0 

202777E 

02 

.9990 

0 

• 235844E 

02 

0 

.237696E 

02 

0 

•239544E 

02 

0 

24 1 385E 

02 

0 

24321 9E 

02 

.9995 

0 

.252604E 

02 

0 

.  25450  7E 

02 

0 

. 256404E 

02 

0 

258294E 

02 

0 

260 1 78E 

02 

.9999 

0 

•290785E 

02 

0 

• 292793E 

02 

0 

•294795E 

02 

0 

296788E 

02 

0 

298775E 

02 

39 


8  of  10 

PERCENTAGE  POINTS  OF  THE  CHI-SQUARE  DISTRIBUTION 


P  V 

7.  1 

7.2 

7.3 

7.4 

7.5 

.0001 

0 

•31 4744E- 

00 

0 • 32988 1 E- 

■00 

0 

. 345373E- 

■00 

0 . 36 1 220E- 

■00 

0 

. 37741 8E- 

-00 

.  0005 

0 

•  505695E 

00 

0.52691  IE 

00 

0 

.54851 9E 

00 

0.570514E 

00 

0 

•592892E 

00 

.0010 

0 

• 622PPHE 

00 

0 . 647089E 

00 

0 

.671 990E 

00 

0.697267E 

00 

0 

. 722975E 

00 

.0050 

0 

•  1 0230  IE 

0  1 

0. 105718E 

01 

0 

•  109175E 

ol 

0. 1 1 267 1 E 

01 

0 

. 1 16206E 

01 

.0100 

0 

• 127809E 

01 

0. 131753E 

0  1 

0 

. 1 35736E 

01 

0. 139757E 

01 

0 

• 143816E 

01 

.0250 

0 

. 1 73730E 

01 

0. 1 7850  9E 

01 

0 

.  1 83324E 

01 

0. 188174E 

01 

0 

. 193059E 

01 

.0500 

0 

.222250E 

01 

0 . 227797E 

0  1 

0 

•  233376E 

01 

0 • 238986E 

01 

0 

• 244626E 

01 

.1000 

0 

.  289762E 

0  1 

0.296239E 

0  1 

0 

. 302743E 

01 

0.309272E 

01 

0 

.315826E 

01 

.2000 

0 

. 389867E 

01 

0.39751 9E 

0  1 

0 

.405 1 90E 

01 

0.41 2878E 

01 

0 

.420584E 

01 

.3000 

0 

.  475643E 

0  1 

0.484166E 

01 

0 

.492700E 

ol 

0.501245E 

01 

0 

•509802E 

01 

.4000 

0 

. 558592E 

0  1 

0 . 567866E 

01 

0 

. 577 1 46E 

01 

0 . 586433E 

01 

0 

. 595725E 

01 

.5000 

0 

• 644562E 

0  1 

0 . 654543E 

0  1 

0 

•  664525E 

01 

0 . 674508E 

01 

0 

•  68449 1 E 

01 

.6000 

0 

.73901 5E 

01 

0.749705E 

01 

0 

• 760389E 

01 

0.771069E 

01 

0 

.781 744E 

01 

.7000 

0 

•849800E 

01 

0.861246E 

0  1 

0 

•87268 1 E 

01 

0.884 1 06E 

01 

0 

• 895520E 

01 

.8000 

0 

•992669E 

01 

0. 1 0050  OE 

02 

0 

.101 730E 

02 

0. 1 02960E 

02 

0 

.1041 87E 

02 

.9000 

0 

• 121526E 

02 

0. 122880E 

02 

0 

. 1 24230E 

02 

0. 125578E 

02 

0 

. 126924E 

02 

.9500 

0 

.  142126E 

02 

0. 1 43577E 

08 

0 

•  1 45025E 

0^ 

0. 1 4647  OE 

02 

0 

. 147912E 

02 

.9750 

0 

.  16J667E 

02 

0. 1 6320 1 E 

02 

0 

. 164732E 

02 

0. 1 66260E 

02 

0 

. 1 67783E 

02 

.9900 

0 

. 1 86388E 

02 

0. 1 880 1 8E 

02 

0 

. 1 89644E 

02 

0. 191265E 

02 

0 

.  1 92882E 

02 

.9950 

0 

•204476E 

02 

0.2061 70E 

02 

0 

. 207859E 

02 

0. 209543E 

02 

0 

.21 1222E 

02 

.9990 

0 

•245047E 

02 

0.246869E 

02 

0 

• 248685E 

02 

0 . 250496E 

02 

0 

. 25230 1 E 

02 

.9995 

0 

•262055E 

02 

0 . 263925E 

02 

0 

'*  265790E 

02 

0 . 267649E 

02 

0 

.  26950  1*E 

02 

.9999 

0 

.300755E 

02 

0 . 302727E 

02 

0 

•  304693E 

02 

0 . 306652E 

02 

0 

.308605E 

02 

P  V 

7.6 

7.7 

7.8 

7.9 

8.0 

.0001 

0 

•  393965E- 

00 

0.41 0859E- 

■00 

0 

•428097E- 

00 

0.445676E- 

-00 

0 

. 463594E- 

•00 

.0005 

0 

•61 5649E 

00 

0 . 638780E 

00 

0 

• 662282E 

oo 

0.686150E 

00 

0 

.71 0379E 

00 

.0010 

0 

• 749049E 

00 

0 . 775504E 

00 

0 

•802335E 

00 

0 . 629537E 

00 

0 

.8571 05E 

00 

.0050 

0 

. 1 1 9779E 

01 

0. 1 23390E 

0  1 

0 

.  1 27037E 

01 

0. 13072  IE 

01 

0 

. 1 34441 E 

01 

.0100 

0 

. 14791  IE 

0  1 

0. 152043E 

0  1 

0 

.  15621  IE 

01 

0. 16041 3E 

01 

0 

. 164650E 

01 

.0250 

0 

. 197977E 

01 

0 • 202928E 

0  1 

0 

.20791  IE 

01 

0.21 2927E 

01 

0 

•217973E 

01 

.0500 

0 

.250296E 

01 

0 • 255996E 

0  1 

0 

.261 724E 

01 

0 • 26748  OE 

01 

0 

•273264E 

01 

•  1  000 

0 

•322A05E 

01 

0 . 329007E 

01 

0 

• 335633E 

01 

0 . 342282E 

01 

0 

.348954E 

01 

.2000 

0 

• 4283Q6E 

01 

0.436045E 

0  1 

0 

.  443800E 

01 

0.451 57 1 E 

01 

0 

.459357E 

01 

.3000 

0 

•518369E 

01 

0.526947E 

0  1 

0 

• 535535E 

01 

0.544134E 

01 

0 

•552742E 

01 

.4000 

0 

• 605022E 

01 

0.614325E 

0  1 

0 

. 623633E 

01 

0 . 632946E 

0  1 

0 

•642264E 

01 

.5000 

0 

• 694474E 

0  1 

0.704458E 

0  1 

0 

•714442E 

ol 

0.724427E 

01 

0 

.7344 12E 

01 

.6000 

0 

.7924 1 5E 

01 

0 .80308 1 E 

0  1 

0 

.81 3742E 

01 

0 . 824400E 

01 

0 

•835053E 

01 

.7000 

0 

. 906925E 

01 

0.91831 9E 

0  1 

0 

•929704E 

01 

0.941 080E 

01 

0 

•952446E 

01 

.8000 

0 

. 1 054 13E 

02 

0. 1 06637E 

02 

0 

•  1 07860E 

02 

0. 1 0908 1 E 

02 

0 

.  1  1 0301 E 

02 

.9000 

0 

•  128267E 

02 

0. 1 29608E 

02 

0 

. 1 30946E 

02 

0. 132282E 

02 

0 

. 133616E 

02 

.9500 

0 

.  149350E 

02 

0. 1 50785E 

02 

0 

. 15221 7E 

02 

0. 1 53647E 

02 

0 

. 1 55073E 

02 

.9750 

0 

. 1 69303E 

02 

0. 1 7081 9E 

02 

0 

. 1 7233 1 E 

02 

0. 173840E 

02 

0 

. 1 75345E 

02 

.9900 

0 

. 1 94494E 

02 

0. 1 961 02E 

02 

0 

. 197706E 

02 

0. 199306E 

02 

0 

.200902E 

02 

.9950 

0 

•212897E 

02 

0.214567E 

02 

0 

.216232E 

02 

0.21 7893E 

02 

0 

.219550E 

02 

.9990 

0 

•254100E 

02 

0 • 255894E 

02 

0 

.  257683E 

02 

0 . 259466E 

02 

0 

.261 245E 

02 

.9995 

0 

•27i 348E 

02 

0.2731 89E 

02 

0 

•275025E 

02 

0 • 276855E 

02 

0 

.  278680E 

02 

.9999 

0 

.31Q552E 

02 

0.31 2492E 

02 

0 

•314426E 

02 

0.316354E 

02 

0 

•318276E 

02 

4o 

9  of  10 

PERCENTAGE  POINTS  OF  THE  CHI-SQUARE  DISTRIBUTION 


P  V 

8.1 

8.2 

8.3 

8.4 

8.5 

.0001 

0 

•  48 i 848E- 

00 

0 • 500435E 

00 

0 

•51 9353E 

00 

0 . 538599E 

00 

0 

.5581 70E 

00 

.0005 

0 

.  734967E 

00 

0.759908E 

00 

0 

.785199E 

00 

0.81 0835E 

00 

0 

. 8368 1 4E 

00 

.0010 

0 

• 885035E 

00 

0.91 3322E 

00 

0 

•  941 962E 

00 

0.970950E 

00 

0 

. 1 00028E 

01 

.0050 

0 

.  1 38196E 

01 

0. 1 4 1986E 

01 

0 

. 1 4581 OE 

01 

0. 149668E 

01 

0 

. 1 53559E 

01 

.0100 

0 

. 1 68920E 

01 

0. 1 73224E 

01 

0 

. 17756  IE 

01 

0. 181930E 

01 

0 

•  186330E 

01 

.0250 

0 

.223OSOE 

01 

0.2281 57E 

0  1 

0 

• 233294E 

01 

0.2384&0E 

01 

0 

. 243655E 

01 

.0500 

O 

•279075E 

01 

0.26491 2E 

01 

0 

. 290775E 

01 

0.296654E 

01 

0 

• 302579E 

01 

.  1  000 

0 

• 355647E 

O  1 

0. J62363E 

01 

0 

• 369099E 

01 

0 • 375856E 

01 

0 

• 382634E 

01 

.2000 

0 

.467 1 59E 

01 

0.474975E 

0  1 

0 

•  482806E 

01 

0 . 49065 1 E 

01 

0 

. 4985 1 OE 

01 

»  3uOU 

0 

.561 360E 

0  1 

0.56998OE 

01 

0 

.  57tse>c:5E 

01 

0 . 56727 1 E 

01 

0 

. 595926E 

01 

.4000 

0 

• 65 1 566E 

01 

0.66091 6E 

0  1 

0 

•670249E 

01 

0 . 679566E 

01 

0 

• 688928E 

01 

.5000 

0 

. 744398E 

01 

0 . 754384E 

0  1 

0 

• 764370E 

01 

0. 774357E 

01 

0 

•764344E 

01 

.6000 

0 

•84570 1 E 

01 

0 . U56346E 

0  1 

0 

•  866987E 

01 

0 . 877623E 

01 

0 

. 888256E 

01 

.7000 

0 

.963603E 

01 

0.9751 51E 

0  1 

0 

. 98649 1 E 

01 

0 . 997822E 

01 

0 

.  100914E 

02 

.6000 

0 

. 1 1 1519E 

02 

0.11 2736E 

02 

0 

. 1 1 3951 E 

02 

0. 1 1 5165E 

02 

0 

. 1 16378E 

02 

.9000 

0 

. 1 34947E 

02 

0. 1 36276E 

02 

0 

. 1 37603E 

02 

0. 138929E 

02 

0 

. 1 4025 1 E 

02 

.9500 

0 

. 1 56497E 

02 

0. 15791 7E 

02 

0 

. 1 59336E 

02 

0. 1 6075 1 E 

02 

0 

.  162  164E 

02 

.9750 

0 

. 1 76646E 

02 

0. 1 78347E 

02 

0 

.  1 79842E 

02 

0. 181 335E 

02 

0 

. 182825E 

02 

.9900 

0 

•202494E 

02 

0 • 204083E 

02 

0 

•205667E 

02 

0 . 207248E 

02 

0 

.  208825E 

02 

.9950 

0 

•22 1 202E 

02 

0.222650E 

02 

0 

.  224494E 

02 

0.2261 34E 

02 

0 

•  227770E 

02 

.9990 

0 

. 263O 1 8E 

02 

0 . 264787E 

02 

0 

.266551 E 

02 

0.26831 OE 

02 

0 

•270065E 

02 

.9995 

0 

.260B00E 

02 

0.26231 5E 

02 

0 

.2841 24E 

02 

0 . 285929E 

02 

0 

•287729E 

02 

.9999 

0 

.3201 93E 

02 

0.3221 04E 

02 

0 

•324009E 

02 

0 . 325909E 

02 

0 

•327803E 

02 

P  V 

8.6 

8.7 

6.8 

8.9 

9.0 

.0001 

0 

•  57SQ84E 

00 

0 . 598277E 

00 

0 

•61 8808E 

00 

0.639653E 

00 

0 

•660809E 

00 

.  0005 

0 

.863I 30E 

00 

0 • 68978 1 E 

00 

0 

.916762E 

00 

0.944069E 

00 

0 

.971 699E 

00 

.0010 

0 

.  1 02995E 

01 

0. 1 05996E 

0  1 

0 

.  1 09030E 

01 

0. 1 12096E 

01 

0 

. 1 151 95E 

01 

.0050 

0 

. 1 57482E 

0  1 

0. 161438E 

01 

0 

.  1 65425E 

01 

0. 169444E 

01 

0 

. 1 73493E 

01 

.0100 

0 

. 1 90762E 

01 

0. 195224E 

0  1 

0 

. 19971 6E 

01 

0  •  204239E 

01 

0 

. 208790E 

01 

.0250 

0 

• 248877E 

01 

0.2541 27E 

0  1 

0 

•  259405E 

01 

0 • 264709E 

01 

0 

•270039E 

01 

.0500 

0 

•3085 1 8E 

01 

0.31 448 1 E 

0  1 

0 

• 320468E 

01 

0 . 326478E 

01 

0 

.33251  IE 

01 

.  1  O00 

0 

. 389432E 

01 

0.396249E 

0  1 

0 

•  403086E 

01 

0 . 409942E 

01 

0 

•41 68 1 6E 

01 

.2000 

0 

• 506362E 

01 

0.51 4269E 

0  1 

0 

.522166E 

01 

0 . 530080E 

01 

0 

.538005E 

01 

.3000 

0 

• 604590E 

0  1 

0.61 3263E 

0  1 

0 

.621 944E 

01 

0 . 630633E 

01 

0 

.639331 E 

01 

.4000 

0 

• 698275E 

01 

0 . 707626E 

01 

0 

•71 698 1 E 

01 

0. 726340E 

01 

0 

.735703E 

01 

.5000 

0 

•  794331 E 

01 

0.60431 9E 

0  1 

0 

•81 4307E 

01 

0 • 624295E 

01 

0 

. 834283E 

01 

.6000 

0 

• 898865E 

u  1 

0.90951 OE 

01 

0 

.920132E 

01 

0.930750E 

01 

0 

.941 364E 

01 

.7000 

0 

.  1 02046E 

02 

0. 1 03177E 

02 

0 

.  1 04306E 

02 

0. 1 05435E 

02 

0 

. 1 06564E 

02 

.6000 

0 

. 1 1 7589E 

02 

0. 1 18799E 

02 

0 

. 1 20008E 

02 

0. 121215E 

02 

0 

• 122421E 

02 

.9000 

0 

. 1 4 1 572E 

02 

0. 1 4289 1 E 

02 

0 

.  1 44208E 

02 

0. 1 45523E 

02 

0 

. 146837E 

02 

.9500 

0 

•  1 63574E 

02 

0. 1 6498 1 E 

02 

0 

. 1 66387E 

02 

0. 1 67789E 

02 

0 

. 1691 90E 

02 

.9750 

0 

. 1 8431  IE 

02 

0. 1 85794E 

02 

0 

.  187275E 

02 

0. 188753E 

02 

0 

. 190226E 

02 

.9900 

0 

.21 0399E 

02 

0.211 969E 

02 

0 

.213536E 

02 

0.2151 OOE 

02 

0 

•21 6660E 

02 

.9950 

0 

.229402E 

02 

0.231031E 

02 

0 

•232655E 

02 

0 . 234276E 

02 

0 

• 235893E 

02 

.9990 

0 

•27 1 8 1 5E 

02 

0 • 27356 1 E 

02 

0 

•  275302E 

02 

0.277039E 

02 

0 

.276772E 

02 

.9995 

0 

• 269524E 

02 

0.291314E 

02 

0 

.293 1 OOE 

02 

0.294881 E 

02 

0 

•296658E 

02 

.9999 

0 

• 329693E 

02 

0.331577E 

02 

0 

•333456E 

04 

0 . 335330E 

02 

0 

.3371 99E 

02 

10  of  10 


percentage  points  of  the  chi-square  DISTRIBUTION 
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Print -out  of  FORTRAN  program  No.  1319 
■by 

H.  A.  Richardson 


This  program  was  used  to  compile  table  1  and  part  of  table  2,  Percentage 
points  of  the  incomplete  gamma -function  ratio ,  for  the  following  P  values 


P  =  ^3.95  (0.05)A 
P  =  A. 5  (0.5)38 
P  =  39  (1.0)50 
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COMPILE  RUN  FORTRAN 

C  PERCENTAGE  POINTS  OF  THE  INCOMPLETE  GAMMA-FUNCTION  RATIO  FOR 

C  THE  FOLLOWING  P-VALUES. 

C  P  =  -.95 ( .05 )4 

C  P  =  4  •  5  (  »  5  )  38 

C  P  =  39 ( 1 .0 )50 

C  JOB  NO.  920603-0010  PROGRAM  NO.  1319 

C  PROGRAMMED  BY  H . A . R I CHARDSON  JULY  13.  1  964 

DIMENSION  A ( 429 . 5 ) .  BC195.12) 

REWIND  1 
REWIND  2 
REWIND  3 

READ  INPUT  TAPE  1.  1.  (<A(I,J),  J  =  1,  5).  I  =  1.  429) 

1  FORMAT (  1  OX .5E 1 4 .6 ) 

N  =  0 

DO  2  I  =  1 .  39 
DO  2  J  =  1  .  11 

K  =  13  -  J 

M  =  5*  I 


L  =  M 

— 

4 

N  =  N 

+ 

1 

IB  = 

1 

DO  2 

I  A 

=  L .  M 

» 

< 

CD 

K  ) 

=  A (N.  IB  ) 

IB  = 

I  B 

+  1 

B  (  1  ,  1 

) 

=  -  0.95 

DO  3 

I 

=  2 .  100 

B  (  I  ,  1 

) 

=  B  (  I  -  1  . 

B  (  1  0  1 

.  1 

)  =  4.5 
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DO  4  I  =  102,  192 

4  B (  I  *  1  )  =  B (  I  -  1,1)  +  0.5 

WRITE  OUTPUT  TAPE  2,  5 

50FORMAT  (  1H1  *  1 3X , 95HPERCENT AGE  POINTS  OF  THE  INCOMPLETE  GAM MA -FUNC T I 
ION  RATIO  FOR  THE  P-VALUES  INDICATED  IN  COLUMN  1 ///l HO , 1  OX ,  1  HP, 6X , 
25H1 . 000,4 X,5H0. 998, 4X,5H0. 990,4 X,5H0. 980, 4X,5HQ. 960,4 X,5HQ. 900, 4X, 
35H0.800,4X,5H0.500,4X,5H0.200,4X,5H0.100,4X, 5H0 .0 1 0 ) 

WRITE  TAPE  3.  ((B(I.J),  J  =  1.  12),  I  =  1,  168) 

WRITE  OUTPUT  TAPE  2,  6,  ((B(I,J).  J  =  1,  12),  I  =  1,  168) 

6  FORMAT  ( 1  HO/1  HO, 4X ,  12F9.3/C5X,  12F9.3)  ) 

WRITE  TAPE  3,  ( ( B ( I , J ) ,  J  =  1,  12),  I  -  170,  192,  2) 

WRITE  OUTPUT  TAPE  2,  7,  <(B(I.J>,  J  =  1,  12),  I  =  170,  192,  2) 

7  FORMAT (5X, 12F9. 3 ) 

END  FILE  2 

REW I ND  1 
REWIND  2 
REWIND  3 
TYPE  ICO 

100  FORMAT (24HEND  OF  JOB,  UNLOAD  TAPES) 

TYPE  101 

101  FORMAT (55H0UTPUT  ON  TAPE  11  MUST  BE  PRINTED  UNDER  FORTRAN  CONTROL) 
STOP  3451 


END 
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Print -out  of  FORTRAN  program  No.  0872 

by 

Elgin  G.  Fry 


This  print -out  is  for  the  main  program ,  which  has  been  titled  Hydrology 
study— a  multipurpose  program  for  selected  cumulative  probability- 
d i s  t r ibut i on  analy s  e s . 
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COMPILE  RUN  FORTRAN 
C  JOB  NO  9206030010  PROG  0872 

C  HYDROLOGY  STUDY  SCS 

C  JOB  NUMBER  920603-0010 

C  PROGRAM  NUMBER  0872 

C  WRITTEN  BY  ELGIN  G.  FRY  DP-SRS-USDA 

C  INPUT  TAPE  IS  10  AND  20 

C  OUTPUT  TAPE  IS  11,  12,.  AND  20 

C  REELS  ON  11  AND  12  ARE  FILLED  IN  SEQUENCE 

C  REEL  ON  20  IS  INPUT  TO  PROGRAM  0911 

C  WORK  TAPE  ON  22 

C 

C  END  OF  JOB  CARD--NINES  PUNCHED  IN  COLUMNS  1-30 

C 

OD I  MENS  I  ON  XI 100, 10) , YI 100, 10) , I YR( 10  0) , XLOGI 10  0, 10) 
1  I  YEAR! 100, 10)  , I  START ( 10) ,  IEND(IO),  DAY(IO),  NI10), 
2SUMX  (10),-  SUM  X  X  (  1 0 )  ,  CVL10),  XM(IO),  GM(IO),  SD(IO) 
3BETAI10),  BETA2 ( 100 ,  GAMMA (10),  SQGAMI10),  U(10,ll) 
4TABLE<  180, 12) »  XLN(IO),  PK(li),  XXNI11),  CHANCE(ll) 
EQUIVALENCE  (X.CAPX) 

REWIND  1 
REWIND  2 
REWIND  3 
REWIND  7 
REWIND  8 
L  =  2 
NUM  =  0 

CHANCE ( 1 )  =  0.0 
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CAPXI  100,10), 
OBS (10) , 

YX( 10) , 

P(  10) ,CZ(  10) , 


CHANCE ( 2 ) 


0.2 


CHANCE ( 3 )  =  1.0 

CHANCE  (  4  ).  =  2.0 

CHANCE ( 5 )  =  4.0 

CHANCE ( 6 )  =  10.0 

CHANCE ( 7 )  =  20.0 

CHANCE ( 8 )  =  50.0 

CHANCE ( 9 )  =  80.0 

CHANCE  (10)  =  90.0' 

CHANCE ( 1 1 )  =  99.0 

D  A  V  (  1 ) 

=  1.0 

DAY! 21 

*  3.0 

DAY( 3) 

=  7.0 

D  A  Y  (  4  1 

=  15.0 

DAY( 5) 

It 

o 

• 

o 

0  A  Y  (  6  ) 

=  60.0 

DAY( 7) 

=  90.0 

DAY( 8) 

=  120.0 

DAY( 9) 

=  183.0 

DAY ( 10) 

=  274.0 

PK (  1  )  = 

■  5. 4-8-3-  6. 

PK ( 2 )  = 

2.878 

PK ( 3  )  = 

2. 326 

P  K  (  4  )  = 

2.054 

PK ( 5  )  = 

1.751 

PK l 6 )  = 

1.282 

P  K  (  7  )  = 

0.  842 

PM8)  = 

0.0 

PK ( 9 )  = 

-0.842 

PK( 10) 

=  -1.282 
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U8 


PK (  1 1 )  =  -2.326 

C  READ  IN  TABLE  OF  U  VALUES  FOR  VARIOUS  LEVELS  OF  PROBABILITY 

READ  TAPE  8,  ( I T ABLE ( I , J ) ,  J=.l,12),  1  =  1,168) 

READ  TAPE  8,  ( { T ABLE ( I , J ) ,  J=l,12),  1=169,180) 

C  READ  IN  PARAMETER  CARD 

20READ  INPUT  TAPE  1,  4,  StATl,  S  T  A  T  2  »•  AREA,  {  I  START  (  I  ) .  IEND(I), 

11=1,10) 

4  F0RMAT(2A5,F10. 3,2013) 

IF  lJSTART(l)-99).  5» 174*174 
C  READ  IN  DATA  CARDS 

5  K  =  I  END ( 1 ) 

DO  6  I  =  1  ,-K 

READ  INPUT  TAPE  1,  8,.  IYRII),  (  X I  I  ,J  ) ,  J=l,10) 

6  WRITE  TAPE  9,  IYRII),  (XII, J),  J=l,10) 

8  FORMAT! I  2 , 7F8. 1 , 3F7. 1 ) 

REWIND  9 

DO  170  KK= 1,10 
DO  9  1=1, K 

9  READ  TAPE  9,  IYRII),  (  X  (  I  ,-J  )  ,  J=l,10) 

REWIND  9 


K 1  =, 

I  START ( KK ) 

K2  = 

I  END! KK ) 

K  3  = 

K2  -  1 

IF  KKt)  2,2,10 

C  CHANGE  STORAGE  LOCATIONS 

10  DO  12  l=Kl,K2 
DO  12  J=l,10 
XL  OG ( I , J )  =  0.0 


YII.,  J)  =  X!  I  ,  J) 
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12  I  YEAR  {  I. »  J  )  =  IYRII) 

C  PLACE  Y  VALUES  ORDERED  HIGH  TO  LOW 

DO  16  J=  1  ,  10 

13  DO  16  I =  K l ,  K  3 

IF  (YII.J)  -  Y(I+1VJ))  14,16,16 

14  HOLD  =  Y( I , J) 

Y  (  I.,  J)  =  Y  (  I  +  1  ,  J  ) 

Y  I  I.+ 1 ,  J  )  =  HOLD 

I  HOLD  =  IYEARI I , J) 

IYEAR(I.J)  =  IYEARI  l  +  l,J) 

IYEARI I  ♦  1  *  J )  = I H  OL  D 
GO  TO  13 

16  CONTINUE 

C  CALCULATE  LOG  TO  BASE  10  OF  Y  VALUES 

DO  18  I =K l »  K2 
DO  18  J e  1  *  10 
IF  (Y( I, J) )  18, IB,  17 

17  XLOG  I  I  ,  J  )  =  LOGXF  I  Y  I  I  ,.  J )  ) 

18  CONTINUE 

C  CALCULATE  CAPIX)  VALUES 

IF  (AREA)  19,19,2  2- 

19  DO  20  I *Kl , K2 
DO  20  J=  1 ,  10 

20  CAPXI  I , J)  =  YI  I i J) 

GO  TO  26 

22  CC  =  0.03719008  /  AREA 
DO  24  I *K  1  , K2 
DO  24  J*=l,10 


24  CAPXI I , J ) 


Y  I  I  ,  J!)  *  DAY!  J)  *  CC 
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C  WRITE  HEADER  LINES 

26  WRITE  OUTPUT  TAPE  L,  28,  STATI,  ST  AT  2 ,  AREA,  KK 

280F0RMAT (  1HL , 1 0 X , 6  I HH I GH E 9 T  MEAN  DISCHARGE  GAMMA  FUNCTION  PARAMETERS 

1  FOR  STATION  ,  2  A5 , 4X  ,,6HARE  A  =  , F 1 0 . 3 , 5X , 3HS E T , I  3 ) 

WRITE  OUTPUT  TAPE  L,  30 

300F0RMAT! 1H0, 115HYR  XII)  LOG  X { 1 )  CAP  X(L)  YR  XI 

12)  LOG  X ( 2 )  CAP  X ( 2 )  YR  X ( 3 )  LOG  X { 3  >  CAP  X  C  3 )  ) 

C  WRITE  DATA  VALUES 

DO  32  I=Ki,K2 

320WRITE  OUTPUT  TAPE  L,  34,  I  Y  E  A  R  {  I  1  )  ,  Y(I,1),  XL0G(I,1),  CAPXII.l), 
1 1  Y  E  AR (  I  ,  2 )  ,  Y ( I i 2 ) »  XL0G(I,2),  C  AP  X (  1,2),  I  YE AR ( I , 3 ) ,  Y(I,3), 

2  XLOG (1,3),  CAPX ( I ,3 ) 

34  FORMAT!  I3,F12.4,F10.6,FL0.4,2(  I8,F12.4,F10.6»F10.4)  ) 

WRITE  OUTPUT  TAPE  L,  36 

360F0RMAT!  1H0, 115HYR'  X(4)  LOG  X  ( 4 )  CAP  X(4)  YR  XI 

15)  LOG  X ( 5 )  CAP  X I  5 )  YR  X(6)  LOG  X(6)  CAP  X I  6 )  ) 

C  WRITE  DATA  VALUES 

DO  38  I =K 1 , K2 

380WRITE  OUTPUT  TAPE  L,  34,  IYEAR(I,4),  Y(I,4),  XL0GII.4),  CAPXII.4), 
1IYEAR!  1,5),  Y(I,5),  XL0G(I,5),  CAPX(I,.5),  IYEAR(I,6),  Y(I,6), 

2XL0G  (  I  ,  6  )  ,  CAPXII.,6) 

WRITE  OUTPUT  TAPE  L,  40 

400F0RMA  T I 1H0»115HYR  X17)  LOG  X(7)  CAP  X(7)  YR  X( 

18)  LOG  X I  8 )  CAP  X  t  8 )  YR  X(9)  LOG  X(9)  CAP  X(9)  ) 

C  WRITE  DATA  VALUES 

DO  42  I  =K l *  K2 

420WRITE  OUTPUT  TAPE  L,  34,  IYEAR(I,7),  Y (  1,7),  XL0G(I,7),  CAPXII.7), 
1 1  YEAR (  1,8),  Y  (  I  ,  8  )  ,  XL0&(I.,8),  CAPXII,8),  IYEAR(I,9),  Y  (  1,9), 


2XL0G  (  I  ,9)  1  CAPXd.,9) 


51 


6  of  15 

WRITE  OUTPUT  TAPE  L»  44 

44  FORMAT! 1H0»35HYR  X I  1 0 )  LOG  X110)  CAP  X (  1 0 )  ) 

DO  46  I  —  K 1 ,  K2 

460WRITE  OUTPUT  TAPE  L ,  34^  IYEAR(I.IO),  Y(I, 10),  XL0G!I,10),  CAPX(I, 
110) 

DO  41  I=K1,K2 
DO  41  J=l,10 
IF  ( C APX ( I »  J ) )  55,55,43 
55  CAPX ( I , J)  =  -9.0 
GO  TO  41 

43  CAPX ( I, J  )  =  LOGXFICAPXl  1, J)  ) 

41  CONTINUE 

KODE  =  K2  -  K1  ♦  1 

WRITE  OUTPUT  TAPE  7,  45,'  ST  ATI  ,  ST  AT  2 ,  AREA,  KK,  KODE 

45  F0RMAT(2A5,.F10.3,  13,  13) 

DO  49  I =K l , K2 

WRITE  OUTPUT  TAPE'  7  ,  47,  (  CAPX  (  I.,.J  )  ,  J=l,.10) 

47  FORMAT! 10F10.4) 

49  CONTINUE 

C  CALCULATE  VALUE  OF  N  FOR  EACH  DAY 

DO  48  1  =  1,  10 

48  Nil)  =  0 

DO  52  I  =  K 1  ,.K2 

DO  52  J=l,10 

IF  (  Y  (  I  ,  J )  )  52,52-,  50 

50  N ( J )  =  N ( J )  ♦  l 
52  CONTINUE 

WRITE  OUTPUT  TAPE  L,  53,  AREA 


53 OF 0RMAT!1HI,30X,66HV0LUME-DURATI0N- PROBABILITY  ANALYSIS  FOR  SELECTE 


52 
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ID  WATERSHEDS  /1H0» 1 4HN AME  OF  S TR E AM , 44 X , 1 3HG A GE  LOCATION/ 

21H0.8HUSGS  NO.  ,  1  5X.,  7HC  TU  NO  .  ,1 1 X  ,.l  6H  PE  R  I  OD  OF  RECORD,  32X, 
3I5HDRAINAGE  AREA  =,F8.3,6H  SQ  MI  / I  HO , 20HDUR AT  I  ON  IN  DAYS  l’.lOX, 
41H3,10X,1H7,9X,2H15,9X,2H30,9X,2H60,9X,2H90,8X,3HI20,8X,3H183,8X, 
53H274  ) 

DO  54  1=1,10 
54  OB  S (  I  )  =  Nil) 

WRITE  OUTPUT  TAPE  L,  56,  (Nil),  1=1,10) 

56  FORMAT ( 10H0N  ,10111) 


CALCULATE  MEAN 

DO  58 

1=1,10 

XM  (  I  ) 

=  0.0 

CVII) 

II 

O 

• 

Q 

CZ  (  I  ) 

II 

O 

• 

Q 

SD  (  I  ) 

=  0.0 

X  L  N  (  I  ) 

=  0.0 

GMI  I  ) 

II 

o 

• 

o 

SUMXI I )  =0.0 

58  SUMXXI I )  =  0.0 
DO  60  I  =  K 1 , K2 
DO  60  J=l,10 

SUMXI  J)  =  SUM  X  (  J ).  +  Y  I  I  ,  J  ) 

60  SUMXXI  J  )  =  SUMXXI  J)  ♦  Y(I,J)  *  Y  (  I  ,.J  ) 

DO  62  J=l,10 

IF  ( OB  S ( J )  )  62,62,61 

61  XM(J)  =  SUM X I J )  /  OB S I J ) 

OCVIJ)  =  (IOBSIJ)  *  OBSIJ)  *  SUMXXI J )  -  OBSIJ)  *  SUMXIJ)  *  SUMXIJ)) 
1/  (IOBSIJ)  -  1.0)  *  SUMXIJ)  *  SUMXIJ)))  **  0.5 


62  CONTINUE 
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C 


WRITE  MEANS,  CV,  AND  SD 


WRITE  OUTPUT  TAPE  L,  64,  (CV(I),  1=1,10) 


64  FORMAT ( 1  OH  CV 


, 1  OF  1 1 . 5 ) 


WRITE  OUTPUT  TAPE  L,  66,  { XM C I ) ,  1  =  1,10) 


66  FORMAT ( 1  OH  X  BAR 


10F11.4) 


C 


CONVERT  LOG  AREA  TO  LOG  OF  X  TO  BASE  E 


DO  70  I =K 1 ,  K2 

DO  70  J*l,10 

IF  ( Y  f I , J ) )  70,70,69 

69  XLOG ( I , J)  =  LOGF ( Y ( I , J  )  ) 

70  CONTINUE 

C  CALCULATE  GEOMETRIC  MEAN  AND  C(Z) 

DO  72  1=1,10 
SUMXXI I )  =  0.0 
72  SUMXI I )  =  0.0 
DO  74  I  =  K  l  ,  K2 
DO  74  J=l,10 

SUMXXIJ)  =  SUMXXIJ)  ♦  XLOG  (  I  ».  J )  *  XLOG(I,J) 

74  SUM X ( J )  =  SUMX(J)  +  XLOGII.J) 

DO  76  J=l,10 
IF  IOBSIJ).)  76,76,75 

750CZIJ)  =  ( (OBSIJ)  *  OBSIJ)  *  SUMXX(J)  -OBSIJ)  *  SUMX(J)  *SUMXIJ)) 
1/  IIOBSIJ)  -  1.0)  *  SUMXIJ)  *  SUMXIJ)))  **  0.5 
SDIJ)  -  SUMXIJ)  /  08 S ( J ) 

XLNIJ)  =  SDIJ) 

CZIJ)  =  CZI  J)  *  ABSFIXLNI J) ) 

GMI J)  =  EXPFI SDI J) ) 

76  CONTINUE 

C  WRI.TE  LOG  GM  AND  THE  GM 
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WRITE  OUTPUT  TAPE  L,  78,  (SD(I),  1=1,10) 

78  F0RMATI10H  LN  G  ,10F11.6) 

WRITE  OUTPUT  TAPE  L,  80,  (GMII),  L=l,10) 

80  FORMATUOH  G  ,10F11.4) 

C  CALCULATE  AND  WRITE  RATIO  OF  X  BAR  /  G 

DO  82  1=1,10 
SD ( I )  =  0.0 
IF  (  OB  S  (  I. )  )  82,82,81 

81  SD (  I  )  =  XM  (  I )  /  GMII) 

82  CONTINUE 

WRITE  OUTPUT  TAPE  L,  84,  ISDII),  1=1,10) 

84  F0RMATI10H  RATIO  X/G,10F11.5) 

C  CALCULATE  AND  WRITE  CV*  AND  Y 

DO  86  1=1,10 

YXI I )  =  0.0 

SUMXl I )  =  0.0 

IF  ( OB  S ( I ) )  86,85,85 

85  SUMXl  I  )  =  (SOU)  *  SDH)  -  1.0)  **  0.5 
YXI I )  =  LOGF ISDII)) 

86  CONTINUE 

WRITE  OUTPUT  TAPE  L,  88,  ISUMXII),.  1=1,10) 

88  FORMAT  I  1  OH  CV*  ,.10F11.5) 

WRITE  OUTPUT  TAPE  L,  90,  IYXII),  1=1,10) 

90  FORMATUOH  Y  ,10F11.6) 

C  CALCULATE  LN  S.D.  AND  LN  S.D.*  (USE  BETA  AREA) 

DO  240  I=.U10 
BETA! I )  =  0.0 
IF  (OBSII))  240,2-40,241 
241  BETA  I  I  )  =  IYXII)  *  2.0)  **  0.5 
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240  CONTINUE 

WRITE  OUTPUT  TAPE  L,  242,  (BETAt I) ,1=1, 10) 

242  FORMAT ( 1  OH  LN  S.D.  ,10F11.6) 

DO  244  1=1,10 

IF  ( OB  S ( I ) )  244,244,243 

243  BETA(I)  =.  BETA(I)  *  (OBSII)  /  (OBS(I)  -  1.0))  **  0.5 

244  CONTINUE 

WRITE  OUTPUT  TAPE  L,  246,  ( BET  A ( I  )  ,  1  =  1,10) 

246  FORMAT { 1  OH  LN  S.O.*  ,10F11.6) 

WRITE  OUTPUT  TAPE  L,  245,  ( C Z (  I  ) ,  1  =  1,10) 

245  FORMATtlOH  LN  S . 0. *  * , 1  OF  1 1 . 6 ) 

C  CALCULATE  BETA  AND  GAMMA  AND  SQR  ROOT  OF  GAMMA 

DO  99  1=1,10 
IF  (OBS(I))  247,247,248 

24 7  GAMMAt I )  =  0.0 
BETAt I )  =  0.0 
SQGAM (  I  )  =  0.0 
GO  TO  99 

248  IF  (YX(I)-0. 5772).  92,94,94 

920GAMM A (  I )  =  (1.0  /  V X ( I  ) >  *  (0.5000876  +  0.  1648852  *  YX (  I  ) 
l-  0.0544274  *  Y X ( I )  *  YX(D) 

GO  TO  98 

94  IF  ( Y  X (  I  )  -  17.0)  96,96,170 

960GAMMA (  I )  =  (8.898919  +  9.059950  *  YX ( I )  +  0.9775373  *  Y X (  I  )  * 

1YX(I))  /  ( YX (  I  )  *  (17.  79728  «■  11.968477  *  YX (  I  )  +  YX (  I )  *  YX I  I  ) )  ) 

98  BETA(I)  =  XM ( I )  /  GAMMA ( I ) 

SQGAM ( I )  =  GAMMA(I)  **  0.5 

99  CONTINUE 


GO  TO  104 
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C  WRITE  GAMMA,  SQR  ROOT  OF  GAMMA  AND  BETA 

104  WRITE  OUTPUT  TAPE  L,  106,  ( GAMMA ( I ) ,  1=1,10) 

106  FORMATtlOH  GAMMA  ,  1  OF  1 1 . 5 ) 

WRITE  OUTPUT  TAPE  L,  108,  (SQGAM(I),  1=1,10) 

108  FORMAT ( 1  OH  SQR  G AMMA , 1  OF  1 1 . 5 ) 

WRITE  OUTPUT  TAPE  L,  110,  (BETA I  I ) ,  1  =  1,10) 

110  FORMATtlOH  BETA  ,10F11.5) 

C  CALCULATE  AND  WRITE  BETA  STAR  (PRIME) 

DO  116  1=1,10 
IF  (AREA)  114,114,112 

112  IF  (  OB  S (  I  )  )  114,114,113 

113  BETA2II)  =  DAY(I).  *  0.03719  *  BET  A  (  I  )  /  AREA 
GO  TO  116 

114  BETA2 ( I )  =  BETA! L) 

116  CONTINUE 

WRITE  OUTPUT  TAPE  L,  118,  ( RET A2 ( I  ) ,  1  =  1,10) 

118  FORMATtlOH  BETA*  ,10F11.5) 

C  SUMX  =  BETA  X  SQR  ROOT  GAMMA 

C  SD  =  BETA*  X  SQR  ROOT  GAMMA 

DO  120  1=1,10 

SUMX ( I )  =  BETA(I)  *  SQGAM(I) 

120  SD (  I )  =  BE  T  A2 ( I )  *  SQGAM(I) 

WRITE  OUTPUT  TAPE  L,  122,  (  SUMX  (  I.)  ,  1=1,10) 

122  FORMATtlOH  B  *  S ( G A ) , 1  OF  1 l . 4 ) 

WRITE  OUTPUT  TAPE  L,  124,  t'SD(l),  1=1,10) 

124  FORMATtlOH  B*  t S ( G A )  , 1  OF  1 1 . 6 ) 

DO  126  1=1,10 
126  P (  I  )  =  GAMMA (  I )  -  l .0 
C  CALCULATE  AND  STORE  U  *  B  *  S  t  G A ) 
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DO  127  1  =  1, .10 
DO  127  J=l,ll 

127  Ul I , J )  =  0.0 
DO  152  J  =  l,.10 

IF  (OBS(J))  152,152,129 

129  IF  C  P  C  J ) -50.0 )  128,128,15 2 

128  DO  1 80  1=2,180 

IF  C  P<  J)-TABLE  (  I  ,.l )  )  131,144,130 

130  CONTINUE 
GO  TO  15? 

131  IF  (  (P(J)-TABLECI-l.l)  )  - l T ABL E (  I , 1 ) -P ( J  )  )  )  132,  132,133 

132  PROB  =  f  P  t  J )  -  T ABLE (1-1,1))  /  ( TABLE!  1,1)  -  TABLE!  1  —  1,11) 

SW  =  0.0 

GO  TO  134 

133  PROB  =  ( TABLE ( I , 1 ) -P ( J ) )  /  f  TABLE (1,1)  -  TABLE  f  1—1,1)) 

SW  =  1.0 

134  DO  142  M  =  1  ,.11 

IF  (TABLE!  1-1, Mi-1)  -  T  ABLE  (  I  ,  M  +  1  )  )  138,  136,  136 

136  IF  (SW)  137,135,137 

137  PROB  T  =  TABLE  (  I  ,  M+!l )  +  (  (T  ABLE  (  I  -  l,Mi  1 )  -  T ABLE (  I , M+ 1 ) )  *  PROB) 

GO  TO  140 

135  PROB  T  =  T  ABLE ( I -1 , M+  1 )  -  (  I  TABLE  I  1-1 ,M+ 1 ) -T ABLE ( I , M  +  l )  )  *  PROB) 
GO  TO  140 

138  IF  (SW)  139,139,141 

139  PROBT=  TABLE! 1-1, M+l)  +  1 ( TABLE < I , M+ 1)  -  T ABLE ( 1-1 , M+l ) )  *  PROB) 
GO  TO  140 

141  PROBT  =  TABLE  (  I  ,  M+l )  -  (I  I  TABLE  (  I  ,M+1  )  -  T  ABLE  (  1-1 ,  M  +  l  )  )  *  PROB) 

140  U ( J  ,M )  =  PROBT*  SD ( J ) 


142  CONTINUE 
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GO  TO  152 

144  DO  146  M=  1 , 1 1 

146  U ( J »  M )  =  T  ABL  E ( I »  M+  1 )  *  SD!J) 

152  CONTINUE 

WRITE  U  *  B  *  S I  3 A )  VALUES 
WRITE  OUTPUT  TAPE  L,  148 

148  FORMAT!  1HO  ,14H  PERCE  NT  CHANCE  ,.43X,.9HU  (  B*S  (  GA  )  ) 

WRITE  OUTPUT  TAPE  L,  154,  (CHANCEII),  (U(J,I),  J=l,10),  1=1,11) 

154  FORMAT  (  1 H  ,.F 9 . 1  , 1  OF  1 1 . 4 ) 

WRITE  OUTPUT  TAPE  L,  161 

1610F0RMAT! 120H0N  =  STATION  YEARS,  CV  =  COEFFICIENT  OF  VARIATION  BY  ME 
1THOD  OF  MOMENTS,  X  BAR  =  ARITHMETIC  MEAN  IN  SECOND-FEET  PER  DAY,  / 
2116H  LN  G  =  NATURAL  LOGARITHUM  OF  THE  GEOMETRIC  MEAN! MEAN  OF  THE  N 
3ATURAL  LOGARITHMS*  OF  THE  ORIGINAL  DATA),  G  =  GEOMETRIC/  120H  MEAN 
4 1 N  SECOND-FEET  PER  DAY,  RATIO  X/G  =  RATIO  OF  ARITHMETIC  MEAN  TO  GE 
50METRIC  MEAN,  CV*  =  COEFFICIENT  OF  VARIATION  BY  /120H  METHOD  OF  MA 
6XIMUM  LIKELIHOOD,  Y  =.  NATURAL  LOGARITHM  OF  THE  RATIO  OF  THE  ARITHM 
7ETIC  MEAN  TO  THE  GEOMETRIC  MEAN,  LN  S.D.  /114H  AND  LN  S.D.*  =  STAN 
8DARD  DEVIATION  OF  THE  NATURAL  LOGARITHMS  FOR  MAXIMUM  LIKELIHOOD,  B 
9IASED  AND  UNBIASED  RESPECT-  ) 

WRITE  OUTPUT  TAPE*  L*  185 

1850F0RMAT  (  120H  IVELY,  GAMMA  =.  SHAPE  STATISTIC,  SQR  GAMMA  =  SQUARE  ROO 
IT  OF  GAMMA,  BETA  A*ND  BETA*  =  RATIO  OF  ARITHMETIC  MEAN  TO  GAMMA  IN/ 
2116H  SECOND-FEET  PER  DAY  AND  VOLUME  IN  INCHES  RESPECTIVELY,  B  *  S( 
3GA )  AND  B*  C  S ( G A )  =  BETA  TIMES  THE  SQUARE  ROOT  OF  /120H  GAMMA 

4 ( STANDARD  DEVIATION)  IN  SECOND-FEET  PER  DAY  AND  VOLUME  IN  INCHES  R 
5ESPECT  I.VELY ,.  U!B*S(GA)  =  VOLUME  IN  INCHES  FOR  /  51H  SELECTED  DUR 
6ATI0NS  AND  PERCENT  CHANCE  AS  INDICATED  /  1H1  ) 


DO  155  1=1,10 


59 


14  of  15 

DO  155  J=1  ,  1 1 
IF  (  U  {  I ,  J  )  )  I  5  1 »  1  5 1  i  1.53 
151  U(  UJ )  =  -9.0 
GO  TO  155 

153  U ( I  ,.  J )  =  LOGXF ( U< I  *  J ).) 

155  CONTINUE 

DO  166  1=1,10 
IF  (OfiSlII)  156,156,157 

157  IF  { GAMMA ( I ) -51 . 0 )  156,156,158 

156  WRITE  OUTPUT  TAPE  7,  164,  I,  (  U  (  I  ,.  J )  ,  J=1,U) 

GO  TO  166 

C  CALCULATE  LOG-NORMAL  SOLUTION 

158  DO  160  J  =  1  , 1 1 

XXN(J)  =  EXPFIXLN(I)  *  PK(J)  *  AftSPI  CZ (  I )  )  ) 

IF  (AREA)  160,160,183 

183  XXN(J)  =  XXN(J)  *  DAY ( I )  *  0.03719  /  AREA 
160  CONTINUE 

WRITE  OUTPUT  TAPE  L,  162,  I,  I  XXNIJ),.  J=  1 ,  11) 

162  FORMAT! 1H0, 12, 4X, 11F10.6) 

DO  163  J=  1  ,  1 1 

IF  ( XXN ( J ) )  171,171,165 

171  XXN ( J )  =  -9.0 
GO  TO  163 

165  XXN(J)  =  LOGXF ( XXN( J) ) 

163  CONTINUE 

WRITE  OUTPUT  TAPE  7,  164,  I,  (XXNIJ),  J=  1 ,  1  1 ) 

164  FORMAT ( I 3,4X, 1 1F10.4) 

166  CONTINUE 

170  NUM  =  NUM  +  (4  *  K2  +.  37) 
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IF  (NUM-20000)  2,172,-172 
172  END  FILE  L 
REWIND  L 
L  =  L  +  1 
NUM  =  0 
GO  TO  2 

174  END  FILE  L 
REWIND  1 
REWIND  L 
KODE  =  999 

WRITE  OUTPUT  TAPE  7,  45,'  ST  AT  1  ,  ST  AT  2 ,  AREA,  KK,  KODE 

END  FILE  7 

REWIND  7 

REWIND  8 

TYPE  176 

176  FORMAT  {  l'OHEND  OF  JOB) 

STOP  9999 


END 
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Print-outs  of  FORTRAN  programs  No.  0910  and  No.  0911 

by 

Elgin  G.  Fry 


In  order  to  plot  on  the  automatic  data  plotter  the  computed  output  from 
program  No.  0872  must  he  coded  and  normal  order  statistics  for  data¬ 
plotter  prepared  (program  No.  0911) . 

Program  No.  0910,  Conversion  of  normal  order  statistics,  was  used  in 
plotting  the  coded  data.  See  the  coded  print-outs  for  Catawba  Creek 
and  Smiths  Fork  forecast,  pages  91 “93  and  98-IOI. 
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compile  run  fortran 

C  CONVERSION  OF  NORMAL  ORDER  STATISTICS 

C  JOS  NUMBER  92u603-O010 

C  PROGRAM  NUMBER  0910 

C  WRITTEN  SY  ELGIN  G.  FRY  DP.  SRS.  USDA 

C 


D I MtNi I  ON 

X  t 

1  00  ) 

» 

A(lo o  ) 

REWIND  1 

REWIND  7 
READ  INPUT 

TAPE 

1 

9 

2. 

(  (X  (  I 

*  J  )  » 

J  =  2  . 

B  )  » 

I 

=  1 

.4  ) 

FORMAT ( 12X 
READ  INPUT 

.  7F9. 

TAPE 

5  ) 

1 

9 

2. 

dXC  I 

.  J  )  . 

J  =  9, 

15). 

I  = 

1.8) 

READ 

I  NPUT 

TAPE 

1 

9 

2. 

(  (  X  (  I 

.  J ) . 

J  =  1  6 

♦  22  ) 

I 

=1.11) 

READ 

I  NPUT 

TAPE 

1 

9 

2. 

d  X  d 

»  J  )  * 

J  =  23 

.29  ) 

I 

=1.15) 

READ 

1  NPUT 

TAPE 

1 

• 

2  * 

d  X  d 

.  J  )  . 

J  =  30 

.36  ) 

I 

=1.18) 

READ 

I  NPUT 

TAPE 

1 

f 

2  . 

d  X  d 

*  J  )  « 

J  =  37 

.43  ) 

I 

= 1 .22 ) 

READ 

I  NPUT 

TAPE 

1 

9 

2  « 

(  (X  (  1 

•  J ) . 

O' 

II 

.50  ) 

I 

=1 .25) 

READ 

I  NPUT 

TAPE 

1 

9 

2  . 

d  X  d 

*  J  )  . 

J  =  5 1 

.57  ) 

I 

= 1 .29 ) 

READ 

I  NPUT 

TAPE 

1 

9 

2  • 

d  X  d 

*  J  )  . 

j  =  sa 

.64  ) 

I 

=1.32) 

READ 

I  NPUT 

TAPE 

1 

9 

2  . 

l  (X  (  I 

.  J). 

J  =  65 

.71  ) 

I 

=1.36) 

REAu 

1  NPUT 

TAPE 

1 

9 

2. 

d  X  (  I 

.  J ) . 

J  =  72 

.78  ) 

I 

=1.39) 

READ 

I  NPUT 

TAPE 

1 

♦ 

2  , 

(  (  X  (  I 

.  J  )  » 

J  =  79 

.85  ) 

I 

=1 .43) 

READ 

I  NPUT 

TAPE 

1 

• 

2  . 

(  (  X  (  I 

*  J  )  . 

J  =  86 

.92  ) 

1 

=1.46) 

READ 

I  NPUT 

TAPE 

1 

9 

2  . 

(  (X  (  I 

.  J) , 

J  =  93 

.99  ) 

I 

=1.50) 

READ 

INPUT 

TAPE 

1 

9 

4  . 

(X  (  I  . 

1  00) 

.  1  =  1 

.50  ) 

4  FORMAT ( 12X.F9. 5 ) 

DO  14  I =2  t 1 00 
DO  6  J  =  1 . I 
6  A ( J  )  -  0.0 

K  =  I 
M  =  I  /  2 
DO  8  J=1  » M 
A  (  J  )  =  X  l  J  .  1  ) 

A (K )  =  -X ( J.  I  ) 

K  =  K  -  1 
8  CONTINUE 
DO  10  J=  1  .  I 
10  A ( J )  =  A ( J  )  +  4.0 

WRITE  OUTPUT  TAPE  7.  12.  (A(J),  J=l.i) 

12  FORMAT ( 1  OF  10. 5) 

14  CONTINUE 
END  FILE  7 
REWIND  1 
REWIND  7 

Type  16 

16  FORMAT ( 1 OHEND  OF  JOS) 

STOP  9999 
END 
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COMPILE  RUN  FORTRAN 

PREPARATION  OF  NORMAL  ORDER  STATISTICS  FOR  DATAPLOTTER 
JOB  NUMBER  920603-0010 
PROGRAM  NUMBER  0911 

WRITTEN  BY  ELGIN  G.  FRY  DP ,  SRS 
DIMENSION  X(5050)  i  A ( 1 00  )  ,  P(ll), 


USDA 

CAP (lOO.lO)i 


SOL (10.11).  N ( 1 0  ) 


REWIND  1 
REWIND  2 
REWIND  3 
REWIND  7 
REWIND  8 
LINE  =  0 
L  =  1 

DO  6  1=2.100 

READ  INPUT  TAPE  7,  2. 

FORMAT (  1  OF  1 0 .5 ) 

1  )  /  2  + 
1  )  /  2 


( A ( J )  ,  J=  1  ,  I  ) 


K  = 

I 

* 

( I  - 

NN  = 

I 

* 

( I  + 

J  = 

0 

DO  4 

M 

=  K 

.NN 

U  = 

J 

+ 

1 

4 

X  (  M  ) 

= 

A 

(  J  ) 

6 

CONT 

INUE 

DO  8 

I 

=  2 

,5050 

8 

X  (  I  ) 

= 

X 

(  I  )  * 

P  (  1  ) 

= 

-&■ 

t  9-040- 

P  (2  ) 

= 

0 

.6878 

P  (  3  ) 

= 

0 

.6326 

P  (  4  ) 

= 

0 

.6054 

P  (  5  ) 

= 

0 

.5751 

P  (6  ) 

= 

0 

.5282 

P  (  7  ) 

= 

0 

.4842 

P  (8  ) 

= 

0 

.4000 

P  (9  ) 

= 

0 

.3  158 

0. 1 


<9.  G4&3-  ,  999?) 

fsloOO 


l  u 
12 

1  4 

1  6 


P(10)  =  O'  27/ 

P ( 1 1 )  =  0. 1674 

READ  INPUT  TAPE  8,  12.  STAT1 

FORMAT (2A5.F10.3.I3.I3) 

IF  (KODE-999)  14.50,50 

DO  18  1=1 , KODE 

READ  INPUT  TAPE  8,  16.  (CAP(I,J), 

FORMAT (  1 0F1 0.4  ) 


STAT2.  AREA,  KK .  KODE 


J=1 ,10) 
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18  CONTINUE 

READ  INPUT  TAPE  8,  20,  ((S0L(I,J),  J= 1 , 1 1 ) ,  1=1,10) 

20  FORMAT (7X , 1 1 F 1 0.4 ) 

DO  22  1=1,10 

22  N (  I  )  =  0 

DO  26  1=1 , KODE 
DO  26  J= 1 , 1 0 

IF  ( CAP ( I , J ) +9 • 0 )  24,26,24 
24  N ( J )  =  N(J)  +  1 

26  CONTINUE 

DO  27  1=1 , KODE 
DO  27  J= 1  ,1  0 

27  CAP (  I  , J  )  =  ( CAP (  I  , J )  +  4.0)  *  0.1 

DO  28  1=1,10 

DO  28  J= 1  ,  1  1 

28  SOL ( I , J )  =  ( SOL ( I , J )  +  4.0)  *  Q»\ 

DO  44  1=1,10 

IF  (NCI))  44,44,30 
30  K  =  N ( I )  *  (NCI)  -  1)  /  2  +  1 

NN  =  N (  I  )  *  ( N (  I  )  +  1)  /  2 

J  =  1 

WRITE  OUTPUT  TAPE  L,  32,  STATU  ST  AT2 ,  KK ,  I,  X(K),  CAP(J,I),  J 
32  FORMAT (2A5»2I3,2F8.4,19x»lH9,I8) 

K  =  K  +  1 
DO  36  M=K,NN 
J  =  J  +  1 

WRITE  OUTPUT  TAPE  L,  34,  STAT1,  STAT2,  KK ,  I,  X(M),  CAP(J,I),  J 
34  FORMAT (2A5,2 I  3 ,2F8. 4 , 1 9x, 1 H8 , 18 ) 

36  CONTINUE 
J  =  J  +  1 

WRITE  OUTPUT  TAPE  L,  38,  STAT1,  STAT2 ,  KK ,  I,  X ( NN  )  ,  CAP(J-1,I),  J 

38  FORMAT ( 2A5 , 2  I  3 , 2F8 . 4 , 1 9x , 1H7, 18 ) 

DO  39  JJ=  1,11 

IF  ( SOL ( I , J J ) )  39,39,41 
41  J  =  J  +  1 
I  I  =  JJ 

WRITE  OUTPUT  TAPE  L,  40,  STAT1,  STAT2,  KK ,  I,  P ( J J ) ,  SOL ( I , J J ) ,  J 
40  FORMAT ( 2A5 , 2  I  3 , 2F8 . 4 , 20x ,  18 ) 

39  CONTINUE 
J  =  J  +  1 

WRITE  OUTPUT  TAPE  L«  38,  STAT1,  STAT2,  KK ,  I,  P(II),  SOL (  I  ,  I  I  )  ,  J 
LINE  =  LINE  +  N ( I )  +  20 
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IF  (LINE-25000)  44.42*42 
42  END  FILE  L 
REWIND  L 
L  =  L  +  1 
LINE  =  0 
44  CONTINUE 
GO  TO  10 
50  REWIND  7 
REWIND  8 
END  FILE  L 
REWIND  L 
TYPE  52 

52  FORMAT ( 1 OHEND  OF  JOB) 
STOP  9999 
END 
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Preparation  of  data  furnished  hy 
U.S.  Geological  Survey  (USGS) 


USGS  provides  the  annual  (water  year)  high-flow  volumes  for  periods  of 
1,  3,  7,  15,  30,  60,  90,  120,  150,  183,  and  2jb  days  for  stream-gage 
records  requested  by  the  Central  Technical  Unit,  SOS. 

The  data  received  from  USGS  are  on  tape  or  ozalid  prints.  If  ozalid 
prints  are  received,  the  IBM  cards  are  key-punched  and  verified  before 
processing.  For  tape,  the  following  programs.  Compile  run  autocoder  by 
P.  Antoinette  Vann  and  format  USGS  data,  FORTRAN  program  No.  1183,  are 
used  to  punch  cards  directly  from  USGS  tape.  The  header  cards  are  in¬ 
serted  before  listing  for  visual  inspection  and  analysis  on  the  com¬ 
puter.  Catawba  Creek  near  Catawba,  Va. ,  Station  No.  02-0185-00,  drain¬ 
age  area  =  3^.0  sq.  mi.,  for  19^-62,  is  a  typical  example. 
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PARflM I  97 1  0 1  OP 

formats  data  from  gs 

JOB  NUMBER  920603-0030 
C  PROGRAM  NUMBER  1183 

D I  MENS  ION  X (  1  O ) 

1  PFAD  INPUT  TAPE  1 ,  6,  K,  J,  N 

6  FORMAT (  I  2 ,  I  8 ,  1 X ,  I  8 ) 

IF  ( K  —  99 )  32,30*30 

30  IF  (J-9999999Q)  A?, 20, 20 
32  IF  ( K— 1  )  34 ,42,42 

34  IF  ( N— 1  )  36,42,42 

36  PRINT  38,  X,  J 
38  FORMAT ( 1 H 1  ,  I  2 ,  I  8 ) 

PUNCH  40,  J 

40  FORMAT (12,18) 

GO  TO  1 

42  BACKSPACE  1 

2  READ  INPUT  TAPE  1,  4,  K,  (X(I),  1=1,10) 

4  FORMAT ( 12, 7F8.2,3F7.2 ) 

PRINT  16,  K ,  ( X (  I  )  ,  1  =  1,10) 

16  FORMAT ( 1  HO ,  I  2 , 7F8. 1 , 3F7. 1  ) 

PUNCH  54,^,  (X(I),  1=1,10) 

54  FORMAT (  12, 7F  8»1,3F7,1) 

GO  TO  1 
20  REWIND  1 
PRINT  22 

22  FORMAT ( 1  HI  ,25X,  1 OHEND  OF  RUN  ////) 

STOP  999 
END 
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Automatic  data  plotting 

The  following  print -out  of  coded  raw  and  computed  data  for  /Catawba  Creek 

and  Smiths  Fork  is  from  tape.  The  automatic  data  plotter^  plots  from  cards 

punched'  from  this  tape  and  a  range  card  made  up  manually. 

1.  Column  1  is  the  station  number  (can  be  any  10  digits  or  symbols). 

2.  Column  2  is  the  set  number  (can  be  10  sets,  only  one  in  this  print-out). 

3.  Column  3  is  the  curve  number  or  field  designation,  which  goes  from  1 
to  10  if  all  fields  are  used.  Both  of  our  examples  contain  10 
columns  of  data. 

4.  Column  4  is  the  coded  plotting  position,  which  is  a  function  of  the 
sample  size  except  as  noted  in  item  6. 

5-  Column  5  gives  the  coded  raw  and  computed  data.  See  item  6. 

6.  Column  6  is  a  code  used  for  line  or  dot  mode  plotting — 9  is  the  code 

on  card  1,  8  on  cards  2  through  19  (sample  size  is  19  years),  and  7 
a  stop  code  on  card  20.  Cards  21  through  31  are  for  the  11  computed 
points  with  their  respective  coded  plotting  positions  given  in  column 
4  (0.9483  is  the  0.0  probability  level  and  0.1674  is  the  99-0  percent 
level) . 

7.  The  last  column  gives  the  card  numbers  for  each  column  or  field.  The 
total  is  the  sample  size  plus  13  (ll  probabilities  and  2  stop  cards 
coded  7)*  Catawba  Creek,  N  =  19  (years  of  record)  plus  13,  has  a 
total  of  32  cards  for  each  of  the  10  curves.  Smiths  Fork,  N  =  l6 
(years  of  record)  plus  13,  has  a  total  of  29  cards  for  all  curves 
except  curve  3  for  which  N  =15,  a  total  of  28  cards. 

The  ADP  plotting  on  log  normal  probability  paper  for  Catawba  Creek 

follows  the  print -out. 
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Dataplotter 


Long  Branch,  N 


Model  5300  manufactured  by  Electronics  Associates, 

J. 


Inc. , 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

L 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

_L 

l 

l 

l 

l 

l 

l 

l 

l 

l 

l 

l 

1 

l 

l 

1 

l 

l 

1 

l 

l 

l 

l 


1 

0.5844 

0.4227 

9 

1 

1 

0.5380 

0.4182 

8 

2 

1 

0.5099 

0.4072 

8 

3 

1 

0.4886 

0.4043 

8 

4 

1 

0.4707 

0.4025 

8 

5 

1 

0.4548 

0.3972 

8 

6 

1 

0.4402 

0.3940 

8 

7 

1 

0.4264 

0.3923 

8 

8 

1 

0.4131 

0.3865 

8 

9 

1 

0.4000 

0.3844 

8 

10 

1 

0.3869 

0.3833 

<L> 

43 

(L) 

rCj 

8 

11 

1 

0.3736 

0.3831 

o 

o 

8 

12 

1 

0.3598 

0.3806 

<D 

•H 

•H 

8 

13 

1 

0.3452 

0.3770 

rH 

8 

14 

1 

0.3293 

0.3704 

u 

O 

ctf 

■g 

o 

8 

15 

1 

0.3114 

0.3613 

|| 

|| 

8 

16 

1 

0.2901 

0.3595 

8 

17 

1 

0.2620 

0.3523 

5 

1 

o 

rH 

01 

o 

1 — 1 

8 

18 

1 

0.2156 

0.3501 

O 

b£) 

£ 

8 

19 

1 

0.2156 

0.3501 

ii 

II 

7 

20 

1 

0.9483 

0.4720 

-p 

-p 

•rH 

a 

21 

1 

0.6878 

0.4350 

-P 

a 

22 

1 

0.6326 

0.4272 

O 

i — 1 

i — 1 

23 

1 

0.6054 

0.4231 

i — 1 

ft 

0 

£ 

Jh 

0 

ft 

24 

1 

0.5751 

0.4185 

u 

25 

1 

0.5282 

0.4108 

tin 

rH 

26 

1 

0.4842 

0.4032 

(D 

bO 

o 

o 

o 

K~\ 

O 

27 

1 

0.4000 

0.3872 

o 

VO 

28 

1 

0.3158 

0.3690 

Note : 

Ph 

29 

1 

0.2718 

0.3584 

30 

1 

0. 1674 

0.3300 

31 

1 

0.1674 

0.3300 

7 

32 

2 

0.5844 

0.4427 

9 

1 

2 

0.5380 

0.4332 

8 

2 

2 

0.5099 

0.4254 

8 

3 

2 

0.4886 

0.4233 

8 

4 

2 

0.4707 

0.4209 

8 

5 

2 

0.4548 

0.4185 

8 

6 

2 

0.4402 

0.4134 

8 

7 

2 

0.4264 

0.4109 

8 

8 

2 

0.4131 

0.4058 

8 

9 

2 

0.4000 

0.4055 

8 

10 

2 

0.3869 

0.4046 

8 

11 

2 

0.3736 

0.4029 

8 

12 

2 

0.3598 

0.4020 

8 

13 

2 

0.3452 

0.4020 

8 

14 

2 

0.3293 

0.4005 

8 

15 

2 

0.3114 

0.3874 

8 

16 

2 

0.2901 

0. 3849 

8 

17 

2 

0.2620 

0.3811 

8 

18 

2 

0.2156 

0.3731 

8 

19 

2 

0.2156 

0.3731 

7 

20 

2 

0.9483 

0.4850 

21 

2 

0.6878 

0.4508 

22 

79 


02-0185.00 

1 

2 

0.6326 

0.4437 

23 

02-0185.00 

1 

2 

0.6054 

0.4401 

24 

02-0185.00 

1 

2 

0.5751 

0.4359 

25 

02-0185.00 

1 

2 

0.5282 

0.4291 

26 

02-0185.00 

1 

2 

0.4842 

0.4224 

27 

02-0185.00 

1 

2 

0.4000 

0.4084 

28 

02-0185.00 

1 

2 

0.3158 

0.3928 

29 

02-0185.00 

1 

2 

0.2718 

0.3839 

30 

02-0185.00 

1 

2 

0.1674 

0.3603 

31 

02-0185.00 

1 

2 

0.1674 

0.3603 

7 

32 

02-0185.00 

1 

3 

0.5844 

0.4488 

9 

1 

02-0185.00 

1 

3 

0.5380 

0.4467 

8 

2 

02-0185.00 

1 

3 

0.5099 

0.4418 

8 

3 

02-0185.00 

1 

3 

0.4886 

0.4393 

8 

4 

02-0185.00 

1 

3 

0.4707 

0.4368 

8 

5 

02-0185.00 

1 

3 

0.4548 

0.4317 

8 

6 

02-0185.00 

1 

3 

0.4402 

0.4306 

8 

7 

02-0185.00 

1 

3 

0.4264 

0.4301 

8 

8 

02-0185.00 

1 

3 

0.4131 

0.4268 

8 

9 

02-0185.00 

1 

3 

0.4000 

0.4266 

8 

10 

02-0185.00 

1 

3 

0.3869 

0.4215 

8 

11 

02-0185.00 

1 

3 

0.3736 

0.4200 

8 

12 

02-0185.00 

1 

3 

0.3598 

0.4181 

8 

13 

02-0185.00 

1 

3 

0.3452 

0.4161 

8 

14 

02-0185.00 

1 

3 

0.3293 

0.4120 

8 

15 

02-0185.00 

1 

3 

0.3114 

0.4057 

8 

16 

02-0185.00 

1 

3 

0.2901 

0.4033 

8 

17 

02-0185.00 

1 

3 

0.2620 

0.3967 

8 

18 

02-0185.00 

1 

3 

0.2156 

0.3888 

8 

19 

02-0185.00 

1 

3 

0.2156 

0.3888 

7 

20 

02-0185.00 

1 

3 

0.9483 

0.4956 

21 

02-0185.00 

1 

3 

0.6878 

0.4632 

22 

02-0185.00 

1 

3 

0.6326 

0.4567 

23 

02-0185.00 

1 

3 

0.6054 

0.4533 

24 

02-0185.00 

1 

3 

0.5751 

0.4494 

25 

02-0185.00 

1 

3 

0.5282 

0.4431 

26 

02-0185.00 

1 

3 

0.4842 

0.4369 

27 

02-0185.00 

1 

3 

0.4000 

0.4242 

28 

02-0185.00 

1 

3 

0.3158 

0.4101 

29 

02-0185.00 

1 

3 

0.2718 

0.4021 

30 

02-0185.00 

1 

3 

0.1674 

0.3812 

31 

02-0185.00 

1 

3 

0.1674 

0.3812 

7 

32 

02-0185.00 

1 

4 

0.5844 

0.4623 

9 

1 

02-0185.00 

1 

4 

0.5380 

0.4602 

8 

2 

02-0185.00 

1 

4 

0.5099 

0.4582 

8 

3 

02-0185.00 

1 

4 

0.4886 

0.4550 

8 

4 

02-0185.00 

1 

4 

0.4707 

0.4514 

8 

5 

02-0185.00 

1 

4 

0.4548 

0.4478 

8 

6 

02-0185.00 

1 

4 

0.4402 

0.4458 

8 

7 

02-0185.00 

1 

4 

0.4264 

0.4456 

8 

8 

02-0185.00 

1 

4 

0.4131 

0.4435 

8 

9 

02-0185.00 

1 

4 

0.4000 

0.4419 

8 

10 

02-0185.00 

1 

4 

0.3869 

0.4400 

8 

11 

02-0185.00 

1 

4 

0.3736 

0.4385 

8 

12 

80 


02-0185.00 

1 

4 

0.3598 

0.4382 

8 

13 

02-0185.00 

1 

4 

0.3452 

0.4319 

8 

14 

02-0185.00 

1 

4 

0.3293 

0.4283 

8 

15 

02-0185.00 

1 

4 

0.3114 

0.4253 

8 

16 

02-0185.00 

1 

4 

0.2901 

0.4232 

8 

17 

02-0185.00 

1 

4 

0.2620 

0.4224 

8 

18 

02-0185.00 

1 

4 

0.2156 

0.4063 

8 

19 

02-0185.00 

1 

4 

0.2156 

0.4063 

7 

20 

02-0185.00 

1 

4 

0.9483 

0.5058 

21 

02-0185.00 

1 

4 

0.6878 

0.4760 

22 

02-0185.00 

1 

4 

0.6326 

0.4700 

23 

02-0185.00 

1 

4 

0.6054 

0.4670 

24 

02-0185.00 

1 

4 

0.5751 

0.4635 

25 

02-0185.00 

1 

4 

0.5282 

0.4578 

26 

02-0185.00 

1 

4 

0.4842 

0.4523 

27 

02-0185.00 

1 

4 

0.4000 

0.4411 

28 

02-0185.00 

1 

4 

0.3158 

0.4287 

29 

02-0185.00 

1 

4 

0.2718 

0.4218 

30 

02-0185.00 

1 

4 

0.1674 

0.4039 

31 

02-0185.00 

1 

4 

0. 1674 

0.4039 

7 

32 

02-0185.00 

1 

5 

0.5844 

0.4838 

9 

1 

02-0185.00 

1 

5 

0.5380 

0.4771 

8 

2 

02-0185.00 

1 

5 

0.5099 

0.4744 

8 

3 

02-0185.00 

1 

5 

0.4886 

0.4731 

8 

4 

02-0185.00 

1 

5 

0.4707 

0.4723 

8 

5 

02-0185.00 

1 

5 

0.4548 

0.4706 

8 

6 

02-0185.00 

1 

5 

0.4402 

0.4704 

8 

7 

02-0185.00 

1 

5 

0.4264 

0.4640 

8 

8 

02-0185.00 

1 

5 

0.4131 

0.4637 

8 

9 

02-0185.00 

1 

5 

0.4000 

0.4584 

8 

10 

02-0185.00 

1 

5 

0.3869 

0.4573 

8 

11 

02-0185.00 

1 

5 

0.3736 

0.4546 

8 

12 

02-0185.00 

1 

5 

0.3598 

0.4501 

8 

13 

02-0185.00 

1 

5 

0.3452 

0.4487 

8 

14 

02-0185.00 

1 

5 

0.3293 

0.4484 

8 

15 

02-0185.00 

1 

5 

0.3114 

0.4454 

8 

16 

02-0185.00 

1 

5 

0.2901 

0 • 4446 

8 

17 

02-0185.00 

1 

5 

0.2620 

0.4352 

8 

18 

02-0185.00 

1 

5 

0.2156 

0.4182 

8 

19 

02-0185.00 

1 

5 

0.2156 

0.4182 

7 

20 

02-0185.00 

1 

5 

0.9483 

0. 5290 

21 

02-0185.00 

1 

5 

0.6878 

0.4973 

22 

02-0185.00 

1 

5 

0.6326 

0.4909 

23 

02-0185.00 

1 

5 

0.6054 

0.4876 

24 

02-0185.00 

1 

5 

0.5751 

0.4838 

25 

02-0185.00 

1 

5 

0.5282 

0.4777 

26 

02-0185.00 

1 

5 

0.4842 

0.4717 

27 

02-0185.00 

1 

5 

0.4000 

0.4593 

28 

02-0185.00 

1 

5 

0.3158 

0.4457 

29 

02-0185.00 

1 

5 

0.2718 

0.4380 

30 

02-0185.00 

1 

5 

0.1674 

0.4179 

31 

02-0185.00 

1 

5 

0. 1674 

0.4179 

7 

32 

02-0185.00 

1 

6 

0.5844 

0.4979 

9 

1 

02-0185.00 

1 

6 

0.5380 

0.4934 

8 

2 

8l 


02-0185.00 

1 

6 

0.5099 

0.4924 

8 

3 

02-0185.00 

1 

6 

0.4886 

0.4907 

8 

4 

02-0185.00 

1 

6 

0.4707 

0.4882 

8 

5 

02-0185.00 

1 

6 

0.4548 

0.4874 

8 

6 

02-0185.00 

1 

6 

0.4402 

0.4870 

8 

7 

02-0185.00 

1 

6 

0.4264 

0.4855 

8 

8 

02-0185.00 

1 

6 

0.4131 

0.4855 

8 

9 

02-0185.00 

1 

6 

0.4000 

0.4787 

8 

10 

02-0185.00 

1 

6 

0.3869 

0.4737 

8 

11 

02-0185.00 

1 

6 

0.3736 

0.4721 

8 

12 

02-0185.00 

1 

6 

0.3598 

0.4721 

8 

13 

02-0185.00 

1 

6 

0.3452 

0.4703 

8 

14 

02-0185.00 

1 

6 

0.3293 

0.4633 

8 

15 

02-0185.00 

1 

6 

0.3114 

0.4591 

8 

16 

02-0185.00 

1 

6 

0.2901 

0.4574 

8 

17 

02-0185.00 

1 

6 

0.2620 

0.4508 

8 

18 

02-0185.00 

1 

6 

0.2156 

0.4392 

8 

19 

02-0185.00 

1 

6 

0.2156 

0.4392 

7 

20 

02-0185.00 

1 

6 

0.9483 

0.5459 

21 

02-0185.00 

1 

6 

0.6878 

0.5143 

'22 

02-0185.00 

1 

6 

0.6326 

0.5080 

23 

02-0185.00 

1 

6 

0.6054 

0.5047 

24 

02-0185.00 

1 

6 

0.5751 

0.5010 

25 

02-0185.00 

1 

6 

0.5282 

0.4950 

26 

02-0185.00 

1 

6 

0.4842 

0.4891 

27 

02-0185.00 

1 

6 

0.4000 

0.4769 

28 

02-0185.00 

1 

6 

0.3158 

0.4635 

29 

02-0185.00 

1 

6 

0.2718 

0.4559 

30 

02-0185.00 

1 

6 

0.1674 

0.4362 

31 

02-0185.00 

1 

6 

0.1674 

0.4362 

7 

32 

02-0185.00 

1 

7 

0.5844 

0.5069 

9 

1 

02-0185.00 

1 

7 

0.5380 

0.5027 

8 

2 

02-0185.00 

1 

7 

0.5099 

0.5019 

8 

3 

02-0185.00 

1 

7 

0.4886 

0.5014 

8 

4 

02-0185.00 

1 

7 

0.4707 

0.5002 

8 

5 

02-0185.00 

1 

7 

0.4548 

0.4968 

8 

6 

02-0185.00 

1 

7 

0.4402 

0.4964 

8 

7 

02-0185.00 

1 

7 

0.4264 

0.4947 

8 

8 

02-0185.00 

1 

7 

0.4131 

0.4936 

8 

9 

02-0185.00 

1 

7 

0.4000 

0.4888 

8 

10 

02-0185.00 

1 

7 

0.3869 

0.4848 

8 

11 

02-0185.00 

1 

7 

0.3736 

0.4824 

8 

12 

02-0185.00 

1 

7 

0.3598 

0.4821 

8 

13 

02-0185.00 

1 

7 

0.3452 

0.4811 

8 

14 

02-0185.00 

1 

7 

0.3293 

0.4795 

8 

15 

02-0185.00 

l 

7 

0.3114 

0.4718 

8 

16 

02-0185.00 

1 

7 

0.2901 

0.4665 

8 

17 

02-0185.00 

1 

7 

0.2620 

0.4576 

8 

18 

02-0185.00 

1 

7 

0.2156 

0.4547 

8 

19 

02-0185.00 

1 

7 

0.2156 

0.4547 

7 

20 

02-0185.00 

1 

7 

0.9483 

0.5536 

21 

02-0185.00 

1 

7 

0.6878 

0.5231 

22 

02-0185.00 

1 

7 

0.6326 

0.5170 

23 

02-0185.00 

1 

7 

0.6054 

0.5139 

24 

82 


02-0185.00 

1 

7 

0.5751 

0.5103 

25 

02-0185.00 

1 

7 

0.5282 

0.5045 

26 

02-0185.00 

1 

7 

0.4842 

0.4989 

27 

02-0185.00 

1 

7 

0.4000 

0.4873 

28 

02-0185.00 

1 

7 

0.3158 

0.4746 

29 

02-0185.00 

1 

7 

0.2718 

0.4675 

30 

02-0185.00 

1 

7 

0.1674 

0.4490 

31 

02-0185.00 

1 

7 

0.1674 

0.4490 

7 

32 

02-0185.00 

1 

8 

0.5844 

0.5135 

9 

1 

02-0185.00 

1 

8 

0.5380 

0.5131 

8 

2 

02-0185.00 

1 

8 

0.5099 

0.5064 

8 

3 

02-0185.00 

1 

8 

0.4886 

0.5061 

8 

4 

02-0185.00 

1 

8 

0.4707 

0.5058 

8 

5 

02-0185.00 

1 

8 

0.4548 

0.5057 

8 

6 

02-0185.00 

1 

8 

0.4402 

0.5023 

8 

7 

02-0185.00 

1 

8 

0.4264 

0.5022 

8 

8 

02-0185.00 

1 

8 

0.4131 

0.5014 

8 

9 

02-0185.00 

1 

8 

0.4000 

0.5013 

8 

10 

02-0185.00 

1 

8 

0.3869 

0.4939 

8 

11 

02-0185.00 

1 

8 

0.3736 

0.4917 

8 

12 

02-0185.00 

1 

8 

0.3598 

0.4916 

8 

13 

02-0185.00 

1 

8 

0.3452 

0.4873 

8 

14 

02-0185.00 

1 

8 

0.3293 

0.4824 

8 

15 

02-0185.00 

1 

8 

0.3114 

0.4785 

8 

16 

02-0185.00 

1 

8 

0.2901 

0.4723 

8 

17 

02-0185.00 

1 

8 

0.2620 

0.4642 

8 

18 

02-0185.00 

1 

8 

0.2156 

0.4588 

8 

19 

02-0185.00 

1 

8 

0.2156 

0.4588 

7 

20 

02-0185.00 

1 

8 

0.9483 

0.5624 

21 

02-0185.00 

1 

8 

0.6878 

0.5313 

22 

02-0185.00 

1 

8 

0.6326 

0.5250 

23 

02-0185.00 

1 

8 

0.6054 

0.5218 

24 

02-0185.00 

1 

8 

0.5751 

0.5181 

25 

02-0185.00 

1 

8 

0.5282 

0.5122 

26 

02-0185.00 

1 

8 

0.4842 

0.5064 

27 

02-0185.00 

l 

8 

0.4000 

0.4945 

28 

02-0185.00 

1 

8 

0.3158 

0.4813 

29 

02-0185.00 

1 

8 

0.2718 

0.4739 

30 

02-0185.00 

1 

8 

0.1674 

0.4547 

31 

02-0185.00 

1 

8 

0.1674 

0.4547 

7 

32 

02-0185.00 

1 

9 

0.5844 

0.5281 

9 

1 

02-0185.00 

1 

9 

0.5380 

0.5220 

8 

2 

02-0185.00 

1 

9 

0.5099 

0.5182 

8 

3 

02-0185.00 

1 

9 

0.4886 

0.5157 

8 

4 

02-0185.00 

1 

9 

0.4707 

0.5153 

8 

5 

02-0185.00 

1 

9 

0.4548 

0.5150 

8 

6 

02^0J.  8  5.00 

1 

9 

0.4402 

0.5138 

8 

7 

02-0185.00 

1 

9 

0.4264 

0.5126 

8 

8 

02-0185.00 

1 

9 

0.4131 

0. 5092 

8 

9 

02-0185.00 

1 

9 

0.4000 

0.5092 

8 

10 

02-0185.00 

1 

9 

0.3869 

0.5068 

8 

1  1 

02-0185.00 

1 

9 

0.3736 

0.5023 

8 

12 

02-0185.00 

l 

9 

0.3598 

0.4984 

8 

13 

02-0185.00 

1 

9 

0.3452 

0.4913 

8 

14 

83 


02-0185.00 

1 

9 

0.3293 

0.4888 

8 

15 

02-0185.00 

1 

9 

0.3114 

0.4873 

8 

16 

02-0185.00 

1 

9 

0.2901 

0.4786 

8 

17 

02-0185.00 

1 

9 

0.2620 

0.4728 

8 

18 

02-0185.00 

1 

9 

0.2156 

0.4630 

8 

19 

02-0185.00 

1 

9 

0.2156 

0.4630 

7 

20 

02-0185.00 

1 

9 

0.9483 

0.5771 

21 

02-0185.00 

1 

9 

0.6878 

0.5439 

22 

02-0185.00 

1 

9 

0.6326 

0.5372 

23 

02-0185.00 

1 

9 

0.6054 

0.5337 

24 

02-0185.00 

1 

9 

0.5751 

0.5297 

25 

02-0185.00 

1 

9 

0.5282 

0.5232 

26 

02-0185.00 

1 

9 

0.4842 

0.5168 

27 

02-0185.00 

1 

9 

0.4000 

0.5036 

28 

02-0185.00 

1 

9 

0.3158 

0.4889 

29 

02-0185.00 

1 

9 

0.2718 

0.4805 

30 

02-0185.00 

1 

9 

0.1674 

0.4585 

31 

02-0185.00 

1 

9 

0.1674 

0.4585 

7 

32 

02-0185.00 

1 

10 

0.5844 

0.5435 

9 

1 

02-0185.00 

1 

10 

0.5380 

0.5293 

8 

2 

02-0185.00 

1 

10 

0.5099 

0.5269 

8 

3 

02-0185.00 

1 

10 

0.4886 

0.5234 

8 

4 

02-0185.00 

1 

10 

0.4707 

0.5230 

8 

5 

02-0185.00 

1 

10 

0.4548 

0.5212 

8 

6 

02-0185.00 

1 

10 

0.4402 

0.5199 

8 

7 

02-0185.00 

1 

10 

0.4264 

0.5175 

8 

8 

02-0185.00 

1 

10 

0.4131 

0.5162 

8 

9 

02-0185.00 

1 

10 

0.4000 

0.5144 

8 

10 

02-0185.00 

1 

10 

0.3869 

0.5124 

8 

11 

02-0185.00 

1 

10 

0.3736 

0.5112 

8 

12 

02-0185.00 

1 

10 

0.3598 

0.5068 

8 

13 

02-0185.00 

1 

10 

0.3452 

0.4994 

8 

14 

02-0185.00 

1 

10 

0.3293 

0.4967 

8 

15 

02-0185.00 

1 

10 

0.3114 

0.4921 

8 

16 

02-0185.00 

1 

10 

0.2901 

0.4833 

8 

17 

02-0185.00 

1 

10 

0.2620 

0.4819 

8 

18 

02-0185.00 

1 

10 

0.2156 

0.4681 

8 

19 

02-0185.00 

1 

10 

0.2156 

0.4681 

7 

20 

02-0185.00 

1 

10 

0.9483 

0. 5883 

21 

02-0185.00 

1 

10 

0.6878 

0.5538 

22 

02-0185.00 

1 

10 

0.6326 

0.5467 

23 

02-0185.00 

1 

10 

0.6054 

0. 5430 

24 

02-0185.00 

1 

10 

0.5751 

0. 5388 

25 

02-0185.00 

1 

10 

0.5282 

0.5320 

26 

02-0185.00 

1 

10 

0.4842 

0.5252 

27 

02-0185.00 

1 

10 

0.4000 

0.5110 

28 

02-0185.00 

1 

10 

0.3158 

0.4952 

29 

02-0185.00 

1 

10 

0.2718 

0.4862 

30 

02-0185.00 

1 

10 

0.1674 

0.4622 

31 

02-0185.00 

1 

10 

0.1674 

0.4622 

7 

32 

SM-FK-FC 

1 

1 

0.5766 

0.5542 

9 

1 

SM-FK-FC 

1 

1 

0.5285 

0.5501 

8 

2 

SM-FK-FC 

1 

1 

0.4990 

0. 5496 

8 

3 

SM-FK-FC 

1 

1 

0.4763 

0.5462 

8 

4 

6  of  12 


a 


V 


84 

SM-FK-FC 

1 

1 

0.4570 

0.5435 

8 

5 

SM-FK-FC 

L 

1 

0.4396 

0.5425 

8 

6 

SM-FK-FC 

1 

1 

0.4234 

0.5422 

8 

7  7  of  12 

SM-FK-FC 

1 

1 

0.4077 

0.5408 

8 

8 

SM-FK-FC 

1 

1 

0.3923 

0.5354 

8 

9 

SM-FK-FC 

1 

1 

0.3766 

0.5342 

8 

10 

SM-FK-FC 

1 

1 

0.3604 

0.5342 

8 

11 

SM-FK-FC 

1 

1 

0.3430 

0.5337 

8 

12 

SM-FK-FC 

1 

1 

0.3237 

0.5326 

8 

13 

SM-FK-FC 

1 

1 

0.3010 

0.5301 

8 

14 

SM-FK-FC 

1 

1 

0.2715 

0.5292 

8 

15 

SM-FK-FC 

1 

1 

0.2234 

0.5225 

8 

16-*- 

SM-FK-FC 

1 

1 

0.2234 

0.5225 

7 

17 

SM-FK-FC 

1 

1 

0.9483 

0.5823 

18 

SM-FK-FC 

1 

1 

0.6878 

0.5613 

19 

SM-FK-FC 

1 

1 

0.6326 

0.5574 

20 

SM-FK-FC 

1 

1 

0.6054 

0.5553 

21 

SM-FK-FC 

1 

1 

0.5751 

0.5531 

22 

SM-FK-FC 

1 

1 

0.5282 

0.5495 

23 

SM-FK-FC 

I 

1 

0.4842 

0.5460 

24 

SM-FK-FC 

1 

1 

0.4000 

0.5391 

25 

SM-FK-FC 

1 

1 

0.3158 

0.5318 

26 

SM-FK-FC 

1 

1 

0.2718 

0.5278 

27 

SM-FK-FC 

1 

1 

0.1674 

0.5179 

28 

SM-FK-FC 

1 

1 

0.1674 

0.5179 

7 

29 

SM-FK-FC 

1 

2 

0.5766 

0.6185 

9 

1 

SM-FK-FC 

1 

2 

0.5285 

0.6182 

8 

2 

SM-FK-FC 

1 

2 

0.4990 

0.6173 

8 

3 

SM-FK-FC 

1 

2 

0.4763 

0.6170 

8 

4 

SM-FK-FC 

1 

2 

0.4570 

0.6152 

8 

5 

SM-FK-FC 

1 

2 

0.4396 

0.6100 

8 

6 

SM-FK-FC 

1 

2 

0.4234 

0.6049 

8 

7 

SM-FK-FC 

1 

2 

0.4077 

0.6045 

8 

8 

SM-FK-FC 

1 

2 

0.3923 

0.6009 

8 

9 

SM-FK-FC 

1 

2 

0.3766 

0.5996 

8 

10 

SM-FK-FC 

1 

2 

0.3604 

0.5978 

8 

11 

SM-FK-FC 

1 

2 

0.3430 

0.5949 

8 

12 

SM-FK-FC 

1 

2 

0.3237 

0. 5919 

8 

13 

SM-FK-FC 

1 

2 

0.3010 

0. 5903 

8 

14 

SM-FK-FC 

1 

2 

0.2715 

0.5892 

8 

15 

SM-FK-FC 

1 

2 

0.2234 

0.5690 

8 

16 

SM-FK-FC 

1 

2 

0.2234 

0.5690 

7 

17 

SM-FK-FC 

1 

2 

0.9483 

0.6638 

18 

SM-FK-FC 

1 

2 

0.6878 

0.6355 

19 

SM-FK-FC 

1 

2 

0.6326 

0. 6299 

20 

SM-FK-FC 

1 

2 

0.6054 

0.6271 

21 

SM-FK-FC 

1 

2 

0.5751 

0.6238 

22 

SM-FK-FC 

1 

2 

0.5282 

0.6186 

23 

SM-FK-FC 

1 

2 

0.4842 

0.6135 

24 

SM-FK-FC 

1 

2 

0.4000 

0.6031 

25 

SM-FK-FC 

1 

2 

0.3158 

0.5918 

26 

SM-FK-FC 

1 

2 

0.2718 

0. 5855 

27 

SM-FK-FC 

1 

2 

0.1674 

0.5694 

28 

SM-FK-FC 

1 

2 

0.1674 

0. 5694 

7 

29 

85 


SM-FK-FC 

1 

3 

0.5736 

0.4690 

9 

1 

SM-FK-FC 

1 

3 

0. 5248 

0.4602 

8 

2 

SM-FK-FC 

1 

3 

0.4948 

0.4591 

8 

3  8  of  12 

SM-FK-FC 

1 

3 

0.4715 

0.4580 

8 

4 

SM-FK-FC 

1 

3 

0 . 4  5  L  6 

0.4568 

8 

5 

SM-FK-FC 

1 

3 

0.4335 

0 . 4568 

8 

6 

SM-FK-FC 

1 

3 

0.4165 

0.4556 

8 

7 

SM-FK-FC 

1 

3 

0.4000 

0.4491 

8 

8 

SM-FK-FC 

1 

3 

0.3835 

0.4380 

8 

9 

SM-FK-FC 

1 

3 

C. 3665 

0.4362 

8 

10 

SM-FK-FC 

1 

3 

0.3484 

0.4342 

8 

1  1 

SM-FK-FC 

1 

3 

0.3285 

0.4301 

8 

12 

SM-FK-FC 

1 

O 

D 

0.3052 

0.4301 

8 

13 

SM-FK-FC 

1 

3 

0.2752 

0.4301 

8 

14 

SM-FK-FC 

1 

3 

0. 2264 

0.4230 

8 

15-*- 

SM-FK-FC 

1 

3 

0. 2264 

0.4230 

7 

16 

SM-FK-FC 

1 

3 

0.9483 

0.5114 

17 

SM-FK-FC 

1 

3 

0.6878 

0.4815 

18 

SM-FK-FC 

1 

3 

0.6326 

0.4755 

19 

SM-FK-FC 

1 

3 

0.6054 

0.4725 

20' 

SM-FK-FC 

1 

3 

0.5751 

0. 4690 

21 

SM-FK-FC 

1 

3 

0.5282 

0.4633 

22 

SM-FK-FC 

l 

3 

0.4842 

0.4578 

2  3 

SM-FK-FC 

1 

3 

0.4000 

0.4465 

24 

SM-FK-FC 

1 

3 

0.3158 

0.4342 

25 

SM-FK-FC 

1 

3 

0.2718 

0.4273 

26 

SM-FK-FC 

1 

3 

0. 1674 

0.4093 

27 

SM-FK-FC 

1 

3 

0. 1674 

0.4093 

7 

28 

SM-FK-FC 

1 

4 

0.5766 

0.5512 

9 

1 

SM-FK-FC 

1 

4 

0.5285 

0.5471 

8 

2 

SM-FK-FC 

1 

4 

0.4990 

0.5470 

8 

3 

SM-FK-FC 

1 

4 

0.4763 

0.5458 

8 

4 

SM-FK-FC 

1 

4 

0.4570 

0.5395 

8 

5 

SM-FK-FC 

1 

4 

0.4396 

0.5354 

8 

6 

SM-FK-FC 

1 

4 

0.4234 

0.5354 

8 

7 

SM-FK-FC 

1 

4 

0.4077 

0.5342 

8 

8 

SM-FK-FC 

1 

4 

0.3923 

0.5295 

8 

9 

SM-FK-FC 

1 

4 

0.3766 

0.5277 

8 

10 

SM-FK-FC 

1 

4 

0.3604 

0.5272 

8 

11 

SM-FK-FC 

1 

4 

0.3430 

0.5263 

8 

12 

SM-FK-FC 

1 

4 

0.3237 

0.5255 

8 

13 

SM-FK-FC 

1 

4 

0.3010 

0.5212 

8 

14 

SM-FK-FC 

1 

4 

0.2715 

0.5201 

8 

15 

SM-FK-FC 

1 

4 

0.2234 

0.5170 

8 

1  6  Mr 

SM-FK-FC 

1 

4 

0.2234 

0.5170 

7 

17 

SM-FK-FC 

1 

4 

0.9483 

0.5848 

18 

SM-FK-FC 

1 

4 

0.6878 

0. 5603 

19 

SM-FK-FC 

1 

4 

0.6326 

0.5556 

20 

SM-FK-FC 

1 

4 

0.6054 

0.5532 

21 

SM-FK-FC 

1 

4 

0.5751 

0.5505 

22 

SM-FK-FC 

1 

4 

0.5282 

0.5462 

23 

SM-FK-FC 

1 

4 

0.4842 

0. 5420 

24 

SM-FK-FC 

1 

4 

0.4000 

0.5335 

25 

SM-FK-FC 

1 

4 

0.3158 

0.5245 

26 

86 


SM-FK-FC 

1 

4 

0.2718 

0.5196 

27 

SM-FK-FC 

1 

4 

0.1674 

0.5070 

28 

SM-FK-FC 

1 

4 

0.1674 

0.5070 

7 

29 

SM-FK-FC 

l 

5 

0.5766 

0.4322 

9 

1 

SM-FK-FC 

1 

5 

0.5285 

0.4079 

8 

2 

SM-FK-FC 

1 

5 

0.4990 

0.4041 

8 

3 

SM-FK-FC 

1 

5 

0.4763 

0.4041 

8 

4 

SM-FK-FC 

1 

5 

0.4570 

0.4041 

8 

5 

SM-FK-FC 

1 

5 

0.4396 

0.4000 

8 

6 

SM-FK-FC 

1 

5 

0.4234 

0.3954 

8 

7 

SM-FK-FC 

1 

5 

0.4077 

0. 3903 

8 

8 

SM-FK-FC 

1 

5 

0.3923 

0.3903 

8 

9 

SM-FK-FC 

1 

5 

0.3766 

0.3845 

8 

10 

SM-FK-FC 

1 

5 

0.3604 

0.3845 

8 

1  1 

SM-FK-FC 

1 

5 

0.3430 

0.3845 

8 

12 

SM-FK-FC 

1 

5 

0.3237 

0.3778 

8 

13 

SM-FK-FC 

1 

5 

0.3010 

0.3778 

8 

14 

SM-FK-FC 

1 

5 

0.2715 

0.3778 

8 

15 

SM-FK-FC 

1 

5 

0.2234 

0.3778 

8 

16 

SM-FK-FC 

1 

5 

0.2234 

0.3778 

7 

17 

SM-FK-FC 

1 

5 

0.9483 

0.4639 

18 

SM-FK-FC 

1 

5 

0.6878 

0.4321 

19 

SM-FK-FC 

1 

5 

0.6326 

0.4257 

20 

SM-FK-FC 

1 

5 

0.6054 

0.4224 

21 

SM-FK-FC 

1 

5 

0.5751 

0.4187 

22 

SM-FK-FC 

1 

5 

0.5282 

0.4126 

23 

SM-FK-FC 

1 

5 

0.4842 

0.4066 

24 

SM-FK-FC 

1 

5 

0.4000 

0.3942 

25 

SM-FK-FC 

1 

5 

0.3158 

0.3806 

26 

SM-FK-FC 

1 

5 

0.2718 

0.3729 

27 

SM-FK-FC 

1 

5 

0.1674 

0.3529 

28 

SM-FK-FC 

1 

5 

0.1674 

0.3529 

7 

29 

SM-FK-FC 

L 

6 

0.5  766 

0.5524 

9 

1 

SM-FK-FC 

1 

6 

0.5285 

0. 5487 

8 

2 

SM-FK-FC 

1 

6 

0.4990 

0.5480 

8 

3 

SM-FK-FC 

1 

6 

0.4763 

0.5470 

8 

4 

SM-FK-FC 

1 

6 

0.4570 

0. 5405 

8 

5 

SM-FK-FC 

1 

6 

0.4396 

0.5393 

8 

6 

SM-FK-FC 

1 

6 

0.4234 

0.5367 

8 

7 

SM-FK-FC 

1 

6 

0.4077 

0.5358 

8 

8 

SM-FK-FC 

1 

6 

0.3923 

0.5308 

8 

9 

SM-FK-FC 

1 

6 

0.3766 

0.5301 

8 

10 

SM-FK-FC 

1 

6 

0.3604 

0.5297 

8 

11 

SM-FK-FC 

1 

6 

0.3430 

0. 5279 

8 

12 

SM-FK-FC 

1 

6 

0.3237 

0.5279 

8 

13 

SM-FK-FC 

1 

6 

0.3010 

0.5241 

8 

14 

SM-FK-FC 

1 

6 

0.2715 

0. 5218 

8 

15 

SM-FK-FC 

1 

6 

0.2234 

0.5190 

8 

16 

SM-FK-FC 

1 

6 

0.2234 

0.5190 

7 

17 

SM-FK-FC 

1 

6 

0.9483 

0.5851 

18 

SM-FK-FC 

1 

6 

0.6878 

0.5613 

19 

SM-FK-FC 

1 

6 

0.6326 

0. 5567 

20 

SM-FK-FC 

1 

6 

0.6054 

0.5544 

21 

SM-FK-FC 

1 

6 

0.5751 

0.5518 

22 

SM-FK-FC 

1 

6 

0.5282 

0.5476 

23 

SM-FK-FC 

1 

6 

0.4842 

0.5435 

24 

SM-FK-FC 

1 

6 

0.4C00 

0. 5354 

25 

SM-FK-FC 

1 

6 

0.3158 

0.5266 

26 

SM-FK-FC 

1 

6 

0.2718 

0.5219 

27 

SM-FK-FC 

1 

6 

0. 1674 

0.5098 

28 

SM-FK-FC 

1 

6 

0.1674 

0.5098 

7 

29 

SM-FK-FC 

1 

7 

0 . 5766 

0.5696 

9 

1 

SM-FK-FC 

1 

7 

0.5285 

0.5662 

8 

2 

SM-FK-FC 

1 

7 

0.4990 

0.5662 

8 

3 

SM-FK-FC 

1 

7 

0.4763 

0 . 5661 

8 

4 

SM-FK-FC 

1 

7 

0.4570 

0.5656 

3 

5 

SM-FK-FC 

1 

7 

0.4396 

0.5653 

8 

6 

SM-FK-FC 

1 

7 

0.4234 

0.5652 

8 

7 

SM-FK-FC 

1 

7 

0.4077 

0.5647 

8 

3 

SM-FK-FC 

1 

7 

0.3923 

0.5640 

8 

9 

SM-FK-FC 

1 

7 

0.3766 

0.5639 

8 

1C 

SM-FK-FC 

1 

7 

0.3604 

0.5636 

8 

1  1 

SM-FK-FC 

1 

7 

0.3430 

0.5622 

8 

12 

SM-FK-FC 

1 

7 

0.3237 

0. 5616 

8 

13 

SM-FK-FC 

L 

7 

0.3010 

0.5605 

8 

14 

SM-FK-FC 

1 

7 

0.2715 

0.5602 

3 

15 

SM-FK-FC 

1 

7 

0.2234 

0.5591 

8 

16 

SM-FK-FC 

l 

7 

0.2234 

0.5591 

7 

17 

SM-FK-FC 

1 

7 

C.9483 

0.5790 

18 

SM-FK-FC 

1 

7 

0.6878 

0.5719 

19 

SM-FK-FC 

L 

7 

0.6326 

0.5703 

20 

SM-FK-FC 

1 

7 

0.6054 

0. 5696 

21 

SM-FK-FC 

1 

7 

0.5751 

0.5688 

22 

SM-FK-FC 

1 

7 

0.5282 

0.5675 

23 

SM-FK-FC 

1 

7 

0.4842 

0.5663 

24 

SM-FK-FC 

1 

7 

0.4000 

0. 5640 

25 

SM-FK-FC 

1 

7 

0.3158 

0.5617 

26 

SM-FK-FC 

1 

7 

0.2718 

0. 5605 

27 

SM-FK-FC 

1 

7 

0. 1674 

0.5577 

28 

SM-FK-FC 

1 

7 

0.1674 

0.5577 

7 

29 

SM-FK-FC 

1 

8 

0.5766 

0.6447 

9 

1 

SM-FK-FC 

1 

8 

0.5285 

0.6407 

8 

2 

SM-FK-FC 

1 

8 

0.4990 

0.6400 

8 

3 

SM-FK-FC 

1 

8 

0.4763 

0.6382 

8 

4 

SM-FK-FC 

1 

8 

0.4570 

0.6377 

8 

5 

SM-FK-FC 

1 

8 

0.4396 

0.6354 

8 

6 

SM-FK-FC 

1 

8 

0.4234 

0.6337 

8 

7 

SM-FK-FC 

1 

8 

0.4077 

0.6328 

8 

8 

SM-FK-FC 

1 

8 

0.3923 

0.6316 

8 

9 

SM-FK-FC 

1 

8 

0.3766 

0.6286 

8 

10 

SM-FK-FC 

1 

8 

0.3604 

0.6281 

8 

11 

SM-FK-FC 

1 

8 

0.3430 

0.6279 

8 

12 

SM-FK-FC 

1 

8 

0.3237 

0.6260 

8 

13 

SM-FK-FC 

1 

8 

0.3010 

0.6255 

8 

14 

SM-FK-FC 

1 

8 

0.2715 

0.6253 

8 

15 

SM-FK-FC 

1 

8 

0.2234 

0.6107 

8 

16 

SM-FK-FC 

1 

8 

0. 2234 

0.6107 

7 

17 

SM-FK-FC 

1 

8 

0.9483 

0.6723 

18 

SM-FK-FC 

1 

8 

0.6878 

0.6525 

19 

SM-FK-FC 

1 

8 

0.6326 

0.6488 

20 

SM-FK-FC 

1 

8 

0.6054 

0.6469 

21 

SM-FK-FC 

1 

8 

0.5751 

0.6448 

22 

SM-FK-FC 

1 

8 

0.5282 

0.6415 

23 

SM-FK-FC 

1 

8 

0.4842 

0.6383 

24 

SM-FK-FC 

1 

8 

0.4000 

0.6319 

25 

SM-FK-FC 

1 

8 

0.3158 

0.6252 

26 

SM-FK-FC 

1 

8 

0.2718 

0.6216 

27 

SM-FK-FC 

1 

8 

0.1674 

0.6125 

28 

SM-FK-FC 

1 

8 

0.1674 

0.6125 

7 

29 

SM-FK-FC 

1 

9 

0.5766 

0.7160 

9 

1 

SM-FK-FC 

1 

9 

0.5285 

0.7158 

8 

2 

SM-FK-FC 

1 

9 

0.4990 

0.7154 

8 

3 

SM-FK-FC 

L 

9 

0.4763 

0.7150 

8 

4 

SM-FK-FC 

1 

9 

0.4570 

0.7141 

8 

5 

SM-FK-FC 

1 

9 

0.4396 

0.7126 

8 

6 

SM-FK-FC 

1 

9 

0.4234 

0.7124 

8 

7 

SM-FK-FC 

1 

9 

0.4077 

0.7120 

8 

8 

SM-FK-FC 

1 

9 

0.3923 

0.7119 

8 

9 

SM-FK-FC 

1 

9 

0.3766 

0.7101 

8 

10 

SM-FK-FC 

1 

9 

0.3604 

0.7098 

8 

1  1 

SM-FK-FC 

1 

9 

0.3430 

0.7097 

8 

12 

SM-FK-FC 

1 

9 

0.3237 

0.7093 

8 

13 

SM-FK-FC 

1 

9 

0.3010 

0.7087 

8 

14 

SM-FK-FC 

1 

9 

0.2715 

0. 7080 

8 

15 

SM-FK-FC 

1 

9 

0.2234 

0.7057 

8 

16 

SM-FK-FC 

1 

9 

0.2234 

0.7057 

7 

17 

SM-FK-FC 

1 

9 

0.9483 

0.7285 

18 

SM-FK-FC 

1 

9 

0.6878 

0.7205 

19 

SM-FK-FC 

1 

9 

0.6326 

0.7188 

20 

SM-FK-FC 

1 

9 

0.6054 

0.7180 

21 

SM-FK-FC 

1 

9 

0.5751 

0.7170 

22 

SM-FK-FC 

1 

9 

0.5282 

0.7156 

23 

SM-FK-FC 

1 

9 

0.4842 

0.7142 

24 

SM-FK-FC 

1 

9 

0.4000 

0.7117 

25 

SM-FK-FC 

1 

9 

0.3158 

0.7091 

26 

SM-FK-FC 

1 

9 

0.2718 

0.7077 

27 

SM-FK-FC 

1 

9 

0.1674 

0.7045 

28 

SM-FK-FC 

1 

9 

0.1674 

0. 7045 

7 

29 

SM-FK-FC 

1 

10 

0.5766 

0.5389 

9 

1 

SM-FK-FC 

1 

10 

0.5285 

0.5382 

8 

2 

SM-FK-FC 

1 

10 

0.4990 

0.5378 

8 

3 

SM-FK-FC 

1 

10 

0.4763 

0.5377 

8 

4 

SM-FK-FC 

1 

10 

0.4570 

0.5356 

8 

5 

SM-FK-FC 

1 

10 

0.4396 

0.5299 

8 

6 

SM-FK-FC 

1 

10 

0.4234 

0.5274 

8 

7 

SM-FK-FC 

1 

10 

0.4077 

0. 5258 

8 

8 

SM-FK-FC 

1 

10 

0.3923 

0. 5228 

8 

9 

SM-FK-FC 

1 

10 

0.3766 

0.5199 

8 

10 

SM-FK-FC 

1 

10 

0.3604 

0.5188 

8 

1  1 

SM-FK-FC 

1 

10 

0.3430 

0.5158 

8 

12 

SM-FK-FC 

1 

10 

0.3237 

0.5155 

8 

13 

SM-FK-FC 

1 

10 

0.3010 

0.5117 

8 

14 
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SM-FK-FC 

1 

10 

0.2715 

0.5093 

8 

15 

SM-FK-FC 

1 

10 

0.2234 

0.4954 

8 

16 

SM-FK-FC 

1 

10 

0.2234 

0.4954 

7 

17 

SM-FK-FC 

1 

10 

0.9483 

0.5813 

18 

SM-FK-FC 

1 

10 

0.6878 

0.5545 

19 

SM-FK-FC 

1 

10 

0.6326 

0.5493 

20 

SM-FK-FC 

1 

10 

0.6054 

0.5466 

21 

SM-FK-FC 

1 

10 

0.5751 

0.5436 

22 

SM-FK-FC 

1 

10 

0.5282 

0.5387 

23 

SM-FK-FC 

1 

10 

0.4842 

0.5339 

24 

SM-FK-FC 

1 

10 

0.4000 

0.5243 

25 

SM-FK-FC 

1 

10 

0.3158 

0.5139 

26 

SM-FK-FC 

1 

10 

0.2718 

0.5082 

27 

SM-FK-FC 

1 

10 

0.1674 

0.4935 

28 

SM-FK-FC 

1 

10 

0.1674 

0.4935 

7 

29 

RUNOFF  VOLUME  IN  INCHES  (DEPTH) 
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REDUCED  VARIATE 
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CTU 


HIGHEST  MEAN  DISCHARGE  GAMMA  FUNCTION  PARAMETERS  FOR  STATION  02-0185.00  AREA  =  3*. 000  SET 
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FORTRAN  program  No.  15^6 

by 

Elgin  G.  Fry 


The  program  is  used  for  plotting  selected  computed  output  of  volume - 
duration  curves  for  dataplotter  on  log-log  paper,  Volume -Duration - 
Probability  Curves  on  an  Annual  Basis. 

At  the  present  time,  the  selected  computed  values  are  punched  on  cards 
processed  to  code  the  computed  values  for  plotting  on  the  automatic  data 
plotter.  Modification  of  program  No.  0872  could  give  this  directly  from 
the  main  program,  thus  eliminating  duplication.  The  coding  is  similar 
to  that  used  in  programs  No.  0910  and  0911. 


o  n 
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C 


compile  run  fortran 

JOB  NUMBER  920603-9045  PROG  1546 

volume-duration  curves  for  dataplotter 

WRITTEN  BY  ELGIN  G.  FRY  DP*  SRS,  USDA 
DIMENSION  P(IO).  SOL (4*10) 

REWIND  1 
REWIND  2 
REWIND  3 
REWIND  8 
LINE  =  0 
L  =  1 

P(10)=  0,4438 
P ( 9 )  =  0.4262 
P ( 8 )  =  0.4079 
P ( 7 )  =  0.3954 

P ( 6 )  =  0.3778 
P ( 5 )  =  0.3477 

P ( 4 )  =  0.3176 
P ( 3 )  =  0.2845 

P ( 2 )  =  0.2477 

P ( 1 )  =  0.2000 

10  READ  INPUT  TAPE  8,  12*  STAT1.  S  T AT  2  *  AREA*  KK.  KODE 

12  FORMAT (2A5 *F 1 0.3  *  I  3*  I  3 ) 

IF  ( KODE  —  999 )  14*50,50 

14  READ  INPUT  TAPE  8*  20.  (  ( SOL (  I  * J  )  *  J= 1.10)*  1=1*4) 

20  FORMAT (4X. 10F7. 4 ) 

DO  30  1=1.4 

DO  30  J=1  ,  1  0 
IF  ( SOL ( I ♦ J) )  30.30,28 

28  SOL(I.J)  =  (LOGXF (SOL ( I « J ) )  +  2*0)  *  0*1 
30  CONTINUE 
J  =  0 

DO  48  1=1,4 

KSW  =  1 
DO  39  JJ= 1  ,  1 0 
IF  (SOL(I.JJ))  39,39,41 
41  J  =  J  +  1 
I  I  =  JJ 

IF  (KSW)  60,60,70 
70  KSW  =  0 

WRITE  OUTPUT  TAPE  L,  72,  STAT1,  STAT2,  KK ,  I,  P(JJ), 
72  FORMAT (2A5,2I3,2F8.4, 1 9x . 1H9, 18 ) 


SOL ( I , JJ ) , 


J 
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GO  TO  39 

60  WRITE  OUTPUT  TAPE  L,  AO,  ST  AT  1  ,  STAT2,  KKt  I  *  P(JJ),  SOL  (  I  , J J  )  ,  J 
4  0  FORMAT ( 2 A5 , 2  I  3 , 2F 8 • 4 , 1 9x , 1 H8 , 18) 

39  CONTINUE 
J  =  J  +  1 

WRITE  OUTPUT  TAPE  L »  38,  ST  AT  1  ,  STAT2,  KK  ,  I«  P(II),  SOL  (  I  ,  I  I  )  »  J 
38  FORMAT (2A5,2I3,2F8.4, 1  9x  , 1H7, 18  ) 

48  CONTINUE 

WRITE  OUTPUT  TAPE  L  «  80 
80  FORMAT (1H  ) 

LINE  =  LINE  +  J  +  1 

IF  (LINE-25000)  44,42,42 

42  END  FILE  L 
REWIND  L 
KKK  =  L  +  9 

TYPE  43,  LINE,  KKK 

43  FORMAT ( 16, 16H  RECORDS  ON  UNIT, 13) 

LINE  =  0 

L  =  L  +  1 

44  GO  TO  10 
50  REWIND  8 

END  FILE  L 
REWIND  L 
KKK  =  L  +  9 
TYPE  43,  LINE,  KKK 
TYPE  52 

52  FORMAT ( 1 OHEND  OF  JOB) 

STOP  9999 
END 


RUNOFF  VOLUME  IN  INCHES  (DEPTH) 
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I  3  7  15  30  60  90  120  183  274 


VOLUME  -  DURATION  -  PROBABILITY  CURVES  ON  AN  ANNUAL  BASIS 


HIGHEST  MEAN  DISCHARGE  GAMMA  FUNCTION  PARAMETERS  FOR  STATION  SM-FK-FC  AREA  =  0.000  SET 
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